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ABSTRACT 

We present a catalogue of 99 candidate clusters and groups of galaxies in the redshift range 0.1 
< Zphot < 1-3 discovered in the Spitzer First Look Survey (FLS). The clusters are selected by their 
Rc - 3.6/im galaxy color-magnitude relation using the cluster red sequence algorithm. Spectroscopic 
redshifts from numerous FLS foUowup projects confirm the photometric redshifts of 29 clusters and 
demonstrate that the Rc - 3.6/im red sequence color provides photometric redshifts with an accuracy 
of Az — 0.04 in the redshift range 0.1 < z < 1.0. Using this cluster sample we compute the 3.6/im, 
4.5/xm, 5.8/im, & 8.0/im cluster luminosity functions (LFs). Similar to previous studies, we find that 
for the bands that trace stellar mass at these redshifts (3.6/im, 4.5/im) the evolution in M* is consistent 
with a passively evolving population of galaxies with a high formation redshift {zf > 1.5). Using the 
3.6/im LF as a proxy for stellar luminosity we remove this component from the MIR (5.8/im & 8.0/im) 
cluster LFs and measure the LF of dusty star formation/AGN in clusters. We find that at z < 0.4 the 
bright end of the cluster 8.0/im LF is well-described by a composite population of quiescent galaxies 
and regular star forming galaxies with a mix consistent with typical cluster blue fractions; however, 
at z > 0.4, regular star forming galaxies are insufhcicnt to account for the excess of 8.0/im galaxies, 
and an additional population of dusty starburst galaxies is required to properly model the 8.0/im LFs. 
Comparison to field studies at similar redshifts shows a strong differential evolution in the field and 
cluster 8.0/im LFs with redshift. At z ~ 0.65 8.0/im-detected galaxies are more abundant in clusters 
compared to the field, but thereafter the number of 8.0/im sources in clusters declines with decreasing 
redshift and by z ~ 0.15, clusters are underdense relative to the field by a factor of ^ 5. The rapid 
differential evolution between the cluster and field LFs is qualitatively consistent with recent field 
galaxy studies that show the star formation rates of galaxies in high density environments are larger 
than those in low density environments at higher redshift. 

Subject headings: infrared: galaxies, galaxies: clusters: photometry — evolution — starburst — funda- 
mental parameters 



1. INTRODUCTION 

Since the compilation of the first large samples of 
galaxy clusters almost 50 years ago (Zwicky 1961; Abell 
1958), clusters have been used as fundamental probes of 
the effect of environment on the evolution of galaxies. 
Over this time, our understanding of this phenomenon 
has grown significantly, and a basic picture of the forma- 
tion and evolution of cluster galaxies between < z < 1 
has emerged. Studies of the stellar populations of cluster 
galaxies via the fundamental plane (e.g., van Dokkum et 
al. 1998; van Dokkum & Stanford 2003; Holden et al. 
2005) and the evolution of the cluster color-magnitude 
relation (e.g., Ellis et al. 1997; Stanford et al. 1998; 
Gladders et al. 1998; Blakeslee et al. 2003; Holden et 
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al. 2004; Mei et al. 2006; Homeier et al. 2006; Tran et 
al. 2007) have shown that the majority of stars in clus- 
ter galaxies are formed at high-redshift (z > 2) and that 
most of the evolution thereafter is the passive aging of 
these stellar populations. Studies of the evolution of the 
near-infrared (NIR) luminosity functions (LFs) of clus- 
ters have shown that not only are the stellar populations 
old, but that the bulk of the stellar mass is already as- 
sembled into massive galaxies at high-redshift (e.g., De 
Propris et al. 1999; Toft et al. 2003; StrazuUo et al. 2006; 
Lin et al. 2006; Muzzin et al. 2007a). Furthermore, 
it appears that the cluster scaling relations seen locally 
(z < 0.1, e.g., Lin et al. 2004; Rines et al. 2004; Lin 
et al. 2003), such as the Halo Occupation Distribution, 
Mass-to-Light ratio, and the galaxy number/luminosity 
density profile are already in place by at least z '--^ 0.5 
(e.g., Muzzin et al. 2007b; Lin et al. 2006). 

These studies suggest a picture where the formation of 
the stars in cluster galaxies, as well as the assembly of the 
galaxies themselves occurs at a higher redshift than has 
yet been studied in detail; and that, other than the pas- 
sive aging of the stellar populations, clusters and cluster 
galaxies have changed relatively little since z ~ 1. This 
picture appears to be a reasonable zeroth-order descrip- 
tion of the evolution of cluster galaxies; however, there 
are still properties of the cluster population which can- 
not be explained within this context. In particular, there 
are significant changes in the morphology (Dressier et al. 
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1997; Postman et al. 2005, Smith et al., 2005), color 
(e.g.. Butcher & Oemler 1984; Rakos & Schombert 1995; 
Smail et al. 1998; EUingson et al. 2001; Margoniner et 
al. 2001; Loh et al. 2008) and star-formation properties 
(e.g., Balogh et al. 1999; Dressier et al. 1999; Poggianti 
et al. 1999; Dressier et al. 2004; Tran et al. 2005a; 
Poggianti et al. 2006, although sec Kodama ct al. 2004) 
of cluster galaxies since 2 ~ 1. The fraction of blue, 
star forming galaxies increases from almost zero at z = 
to as much as 50% at z ~ 0.5 (the so-called Butcher- 
Oemler Effect), and correspondingly, the fraction of SO 
galaxies in clusters drops by a factor of 2-3, with similar 
increase in the number of spiral/irregular galaxies over 
the same redshift range (Dressier et al. 1997). Naively, 
these results suggest that gas-rich, star-forming galax- 
ies at high-redshift have their star-formation truncated 
by the cluster environment at moderate redshift and be- 
come the dominate SO population seen locally. How such 
a transformation occurs, and how it avoids leaving a no- 
table imprint on the stellar populations, is still not well- 
understood. 

Citing an abundance of post-starburst (k+a) galaxies 
in clusters at z ~ 0.4, Poggianti et al. (1999) and Dressier 
et al. (2004) suggested that there may be an abundance 
of dusty starburst galaxies in clusters at moderate red- 
shift, and that the dusty starburst and k-|-a galaxies may 
represent the intermediate stages between regular star 
forming late- type galaxies and SO galaxies (e.g., Shioya 
et al. 2004; Bekki & Couch 2003). In particular, they 
suggested that the cluster e(a)^ galaxies would be the 
best candidates for dusty starburst galaxies because their 
inferred star formation rates appear larger from Ha emis- 
sion than from [Oil] emission. If the cluster environment 
excites a dusty starburst from harassment, tidal interac- 
tion, or ram-pressure stripping, then this may quickly 
deplete a star forming galaxy of its gas, transforming it 
first into a k+a galaxy, and then leaving it an SO. More 
detailed work on two z ^ 0.5 clusters by Moran et al. 
(2005) also showed an abundance of starbursting galax- 
ies conspicuously near the cluster virial radius, suggest- 
ing a environmental origin to their "rejuvenation" . ISO 
observations of relatively nearby clusters have detected 
significant amomits of dust-obscured star formation (e.g., 
Fadda ct al. 2000; Due et al. 2002; Biviano et al. 2004; 
Coia et al. 2005), and this has recently been confirmed 
at even higher redshift {z = 0.2 - 0.8) by Spitzer obser- 
vations (Geach et al. 2006; Marcillac et al. 2007; Bai 
et al. 2007; Fadda et al. 2007; Saintongc et al. 2008; 
Dressier ct al. 2008). Despite this, it is currently unclear 
whether there is a population of dusty starbursts which 
is sufficiently abundant to be the progenitors of the large 
number of cluster k-l-a galaxies. 

Alternatively, there is evidence from other cluster sam- 
ples that the SO population may simply be the result 
of the truncation of star formation in infalling late-type 
galaxies via gas strangulation (e.g., Abraham et al. 1996; 
Balogh et al. 1999; Treu et al. 2003; Moran et al. 
2006) and that no accompanying starburst occurs. Most 
likely, the star formation and morphology of galaxies are 

^ An c(a) galaxy is defined as a galaxy with EW([On]) < -SA 
and EW(H5) > dA by Dressier et al. (1999). These are emission 
line galaxies with a strong A star component to their spectrum 
suggesting a recent, possibly obscured, burst of star formation. 



transformed both "actively" (as in a starburst triggered 
from merging/harassment/tidal forces) and "passively" 
(from gas strangulation or ram-pressure stripping), and 
that the magnitude of each effect varies significantly from 
cluster-to-cluster and possibly by epoch, which may ex- 
plain why studies of small numbers of clusters have found 
discrepant results. Interestingly, both processes can be 
active within massive clusters as was demonstrated by 
Cortese et al. (2007), who found two interesting galax- 
ies in Abell 1689 and Abell 2667, one of which seems to 
be undergoing gas strangulation and ram-pressure strip- 
ping, while the other is experiencing an induced star- 
burst. There is evidence that galaxies in clusters that 
are less dynamically relaxed have larger star formation 
rates (e.g., Owen et al. 1999; Metevier et al. 2000; Moss 
& Whittle 2000; Owen et al. 2005; Moran et al. 2005; 
Coia et al. 2005, and numerous others) and that the 
accretion of large substructures induces starbursts from 
harassment and tidal forces. 

The most obvious way to understand whether dusty 
starbursts are important in the evolution of cluster galax- 
ies is to observe their abundances directly in the mid- 
infrared (MIR). In particular, differences in the MIR LFs 
of the cluster and field environments can be used to deter- 
mine if dusty starbursts are more common in the cluster 
environment. If so, it would suggest that environmental 
processes may be responsible for triggering these events. 

The InfraRed Array Camera (IRAC) onboard Spitzer 
provides a unique tool for studying this problem. IRAC 
images in 4 bands simultaneously (3.6/im, 4.5/um, 5.8/xm, 
8.0^m) and this is particularly advantageous because 
3.6/im and 4.5/im observations are a good proxy for the 
stellar mass of cluster galaxies between < z < 1, and 
5.8/im and S.Oiim are sensitive to emission from warm 
dust (i.e., from dusty star forming regions) over the same 
redshift range. In particular, the Polycyclic Aromatic 
Hydrocarbons (PAHs) emit strong line emission at rest 
frame 3.3/^m, 6.2/xm, 7.7/im, 8.6/im, and 11.3/um (e.g., 
Gillett et al. 1973; Willner et al. 1977). These features, 
in addition to the warm dust continuum, are sensitive 
indicators of dusty star formation, and several studies 
have already shown a good correlation between 8.0/xm 
flux and star formation rate (SFR^"; e.g., Calzetti et al. 
2005; Wu et al. 2005; Calzetti et al. 2007). Therefore, 
examining the suite of IRAC cluster LFs at redshifts 
< z < 1 shows both the evolution of the majority of stel- 
lar mass in cluster galaxies, as well as the evolution of 
dusty star formation in the same galaxies. 

The obvious approach to measuring the presence of 
dusty star formation in clusters is to observe a hand- 
ful of "canonical" galaxy clusters with IRAC. However, 
given that determining the LF from a single cluster suf- 
fers significantly from Poisson noise, and perhaps most 
importantly, is not necessarily representative of the av- 
erage cluster population at a given mass/epoch, a better 
approach would be to stack large numbers of clusters 

Although there is a direct correlation between 8fj,m flux and 
SFR, the scatter in the correlation is approximately a factor of 2 for 
metal rich galaxies in the local universe, and metal poor galaxies 
can deviate by as much as a factor of 50 (e.g., Calzetti et al. 2007). 
Because of the large scatter and nietallicity dependence, throught 
this paper we do not use the 8/im data to quantitatively measure 
SFRs. Instead, we use the presence of enhanced 8/im flux as a 
qualitative indicator of increased dusty star formation. 
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in order to improve the statistical errors, and avoid pe- 
culiarities associated with individual clusters. This ap- 
proach requires targeted observations of numerous clus- 
ters, which is time-consuming compared to other alter- 
natives. For example, large-area Spitzer surveys such 
as the 50 deg^ Spitzer Wide-area Infrared Extragalactic 
Survey (SWIRE", Lonsdale et al. 2003), the 8.5 deg^ 
IRAC Shallow Survey (Eisenhardt et al. 2004), and the 
3.8 deg^ Spitzer First Look Survey (FLS-^^, Lacy et al. 
2005) are now, or soon-to-be, publically available and 
these fields already contain significant amounts of opti- 
cal photometry. These wide optical-IRAC datasets can 
be employed to find clusters in the survey area itself using 
optical cluster detection methods such as the cluster red 
sequence (CRS) technique (Gladders & Yee 2000, here- 
after GYOO), or photometric redshifts (e.g., Eisenhardt et 
al. 2008; Brodwin et al. 2006). Subsequently, the IRAC 
survey data can be used to study the LFs of clusters at a 
much larger range of masses and redshifts than could be 
reasonably followed up by Spitzer. Furthermore, these 
surveys also provide panoramic imaging of clusters out 
to many virial radii, something that has thus far rarely 
been attempted because it is time-consuming. 

Finding clusters with the CRS algorithm is relatively 
straightforward with the ancillary data available from 
these surveys. The technique exploits the fact that the 
cluster population is dominated by early type galaxies, 
and that these galaxies form a tight red sequence in color- 
magnitude space. If two filters which span the 4000A 
break are used to construct color-magnitude diagrams, 
early types are always the brightest, reddest galaxies at 
any redshift (e.g., GYOO) and therefore provide signifi- 
cant contrast from the field. The CRS technique is well- 
tested and provides photometric redshifts accurate to ^ 
5% (Gilbank et al. 2007a; Blindert et al. 2004) as weU as 
a low false-positive rate (<5%, e.g., Gilbank et al. 2007a; 
Blindert et al. 2004; Gladders & Yee 2005). The method 
has been used for the 100 deg^ red sequence Cluster Sur- 
vey (RCS-1, Gladders & Yee 2005) and is also being used 
for the next generation, 1000 deg^ RCS-2 survey (Yee et 
al. 2007). Variations of the red sequence method have 
also been used to detect clusters in the Sloan Digital Sky 
Survey (the "BCGmax" algorithm, Koester et al. 2007; 
Bahcall et al. 2003) as well as in the fields of X-ray sur- 
veys (e.g. Gilbank et al. 2004; Barkhouse et al. 2006). 

In this paper we combine the Spitzer FLS Rc-band 
and 3.6/im photometry and use it to detect clusters with 
the CRS algorithm. Given the depth of the data, and 
that the Rc - 3.6/im filter combination spans the rest- 
frame 4000 A break to z > 1, we are capable of detecting 
a richness-limited sample of clusters out to z 1. Using 
the sample of clusters discovered in the FLS we compute 
the 3.6/im, 4.5/im, 5.8/im, and 8.0/im LFs of clusters 0.1 
< z < 1.0 and study the role of dusty star formation in 
cluster galaxy evolution. A second paper on the abun- 
dance of dusty starburst galaxies detected at 24/im in 
the same clusters using the FLS MIPS data is currently 
in preparation by Muzzin et al. (2008). 

The structure of this paper is as follows. In §2 we give 

SWIRE data are publically available at 
|http: //swire.ipac. caltech.ed u/swire/l 

The "FL S dat a are publically available at 
|http://ssc. spitzer . caltech . edu/fls/ 1 



a brief overview of the optical, IRAC, and spectroscopic 
data used in the paper. Section 3 describes the cluster- 
finding algorithm used to detect clusters and §4 contains 
the FLS cluster catalogue, and a basic description of its 
properties. In §5 we present the IRAC cluster LFs and 
§6 contains a discussion of these results as well as a com- 
parison of the cluster and field LFs. We conclude with 
a summary in §7. Throughout this paper we assume an 
— 0.3, flA — 0.7, Ho = 70 km s~^ Mpc ~^ cosmology. 
All magnitudes are on the Vega system. 

2. DATA SET 

2.1. Spitzer IRAC Data and Photometry 

The IRAC imaging data for this project was observed 
as part of the publically available, Spitzer First Look 
Survey (FLS; see Lacy et al. 2005 for details of the 
data acquisition and reduction). The FLS was the first 
science survey program undertaken after the telescope's 
in-orbit-checkout was completed. It covers 3.8 square 
degrees and has imaging in the four IRAC bandpasses 
(3.6/im, 4.5/im, 5.8/im, 8.0/im). The FLS is a shallow 
survey with a total integration time of only 60 seconds 
per pixel. Because IRAC images all four channels si- 
multaneously, the total integration time is identical in 
each channel. The resulting 5cr limiting flux densities are 
20, 25, 100, and 100 /iJy in the 3.6/im, 4.5/im, 5.8/im, 
8.0/im bandpasses, respectively. These flux densities cor- 
respond to Vega magnitudes of 18.0, 17.2, 15.2, and 14.6 
mag, respectively. The 50% completeness limits for the 4 
channels are 18.5, 18.0, 16.0, 15.4 mag and hereafter we 
use these limits for the cluster finding algorithm (§3) and 
computing the cluster LFs (§5). The data was corrected 
for completeness using a third-order polynomial fit to 
the survey completeness as a function of magnitude de- 
termined by Lacy et al. (2005). Lacy et al. compared 
their completeness estimates, made using artifical galax- 
ies, to completeness estimates determined by comparing 
the recovery of sources in the FLS to a deeper "verifi- 
cation strip" . The completeness was similar using both 
methods; however, in some cases the latter suggested it 
might be higher by ^ 10-15%. When counting galaxies 
we have multiplied the formal uncertainties by an addi- 
tional ± 20% of the completeness correction to account 
for this additional uncertainty. 

Photometry for the IRAC data was performed using 
the SExtractor (Berlin & Arnouts 1996) package. For 
each channel, four aperture magnitudes plus an isopho- 
tal magnitude are computed. The four apertures used 
are 3, 5, 10, and 20 IRAC pixels in diameter (3.66", 
6.10", 12.20", 24.40"). The aperture magnitudes are 
corrected for the flux lost outside the aperture due to 
the large diffraction-limit of the telescope and the sig- 
nificant wings of the IRAC point spread function (PSF). 
The aperture corrections are computed from bright stars 
within the FLS field and are listed and discussed fur- 
ther in Lacy et al. (2005). The majority of galaxies 
with 3.6/im > 15.0 mag are unresolved, or only slightly 
resolved at the resolution of the 3.6/im bandpass and 
therefore the 3 pixel aperture corrected magnitude pro- 
vides the best total magnitude. For galaxies which are 
extended and resolved, this small aperture is an under- 
estimate of their total magnitude. For these galaxies, 
a "best" total magnitude is measured by estimating an 
optimum photometric aperture using the isophotal mag- 
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nitudes. The geometric mean radius of the isophote {rm 
= (A/tt)'^ ''^, where A is the isophotal area) is compared 
to the radius of each of the 4 apertures used for the aper- 
ture magnitudes (ri,r2,r3,r4). If < 1.1 rap, then that 
aperture magnitude is chosen as the best total magni- 
tude. For objects with > 1.1 r4 the isophotal mag- 
nitude is used as the best total magnitude. When mea- 
suring the Rc - 3.6y^m colors, we always use the 3 pixel 
aperture-corrected magnitude, even for resolved galaxies 
(see discussion in §2.3). 

Object detection was performed separately in all 4 
channels and these catalogues were later merged using 
a 1.8" search radius. Tests of this matching (Lacy et al. 
2005) show that this radius provides the most reliably 
matched catalogues. 

2.2. Optical Data 

The ground-based Cousins Rc-band (hereafter "R- 
band") imaging used in this study was obtained as part 
of the FLS campaign and is also publically available. R- 
band imaging covering the entire FLS IRAC and MIPS 
fields was observed on the Kitt Peak 4m Mayall tele- 
scope using the MOSAIC-1 camera. MOSAIC-1 consists 
of eight 4096 x 2048 CCDs, and has a field-of-view of 36 
X 36 arcmin with a pixel scale of 0.258 arcseconds per 
pixel. Data reduction was performed using the NOAO 
IRAF mscred package and procedures, and galaxy pho- 
tometry was performed using the SExtractor (Bertin & 
Arnouts 1996) package. Typical seeing for the images 
was ~ 1.1 arcseconds and the 5cr limiting magnitude in 
an aperture of 3 arcseconds is 24.7 mag. The 50% com- 
pleteness limit in the same aperture is ^ 24.5 mag. A 
complete discussion of the data reduction, object finding, 
and photometry can be found in Fadda et al. (2004). 
For this study we performed additional photometry to 
that publically available in order to measure fluxes in 
a slightly larger 3.66" aperture which matches with the 
smallest aperture of the IRAC data (D. Fadda, private 
communication) . 

The mean absolute positional error in the astrometry 
for the R-band data is 0.35" (Fadda et al. 2004) and 
the mean positional error in the astrometry of bright 
(faint) sources in the IRAC catalogue is 0.25" (1.0") 
(Lacy et al. 2005). Given these uncertainties, as well 
as the large IRAC pixel scale, the R-band catalogue was 
matched to the IRAC catalogue by looking for the closest 
object within 1.5 IRAC pixels (1.8") of each IRAC de- 
tection. Tests of matching radii ranging between 0.3 and 
3.0 IRAC pixels (0.37" - 3.66") showed that the number 
of matches increased rapidly using progressively larger 
radii up to ~ 1.5 IRAC pixels and thereafter the gain in 
the number of matches with increasing radius was rela- 
tively modest, suggesting that the majority of additional 
matches were likely to be chance associations. Given 
that the IRAC astrometry is calibrated using bright stars 
from the Two Micron All Sky Survey (2MASS, Strutskie 
ct al. 2006), whereas the R-band data was astrometri- 
cally calibrated using the USNO-A2.0 catalogue (Monet 
et al. 1998) an additional concern was the possibility 
of a systematic linear offset between the two astromet- 
ric systems. We attempted to iteratively correct for any 
systematic offset by shifting the IRAC astrometry by the 
median offset of all matched sources and then rematch- 
ing the catalogues; however, multiple iterations could not 



converge to a solution significantly better than the ini- 
tial 0.2" offset seen between the two systems. Given that 
this offset is less than the quoted positional errors in the 
two systems, it suggests that any systematic offset be- 
tween the 2MASS and USNO-A2.0 system in the FLS 
field is less than the random positional error in the R- 
band and IRAC data themselves. The iterative refine- 
ments increased/decreased the total number of matches 
by ± 0.05 - 0.3% depending on which iteration. Given 
these small variations, and the lack of further evidence 
for a systematic offset between the coordinate systems, 
the final matched catalogue uses the original IRAC and 
R-band astrometry. 

In approximately 4% of cases more than one R-band 
object was located within the search radius. In these 
cases, the object closest to the IRAC centroid was taken 
as the match. The space density of R-band soiirccs is 
approximately 5 times higher than the number of IRAC 
sources at these respective depths. This suggests that at 
most, 20% of R-band sources have an IRAC counterpart 
at the respective depths. 

When there arc multiple R-band matches for an IRAC 
detection, the majority of cases will be where only one of 
the R-band detections is the counterpart of the IRAC 
detection, and our approach will provide correct col- 
ors. Nevertheless, a certain percentage of the multi- 
ple matches will be when two R-band objects, both of 
which have IRAC counterparts, have these counterparts 
blended together into a single IRAC detection due to the 
large IRAC PSF. Because the IRAC source is a blend 
of two objects, but we use only one R-band counterpart, 
these objects will be cataloged as brighter and redder 
than they truly are. However, because only 4% of IRAC 
sources have multiple R-band matches, and the proba- 
bility that both of those R-band sources have an IRAC 
counterpart is roughly, 20%^ = 4%, this suggests that 
only 4% X 4% = 0.16% of all IRAC sources are blended 
sources where only one R-band galaxy has been identi- 
fied as the counterpart. 

Although this estimated contamination is small, clus- 
ters have greater surface densities of galaxies than the 
field, and therefore it might be expected that cluster 
galaxies are blended more frequently than field galax- 
ies. We measiircd the frequency of multiple matches for 
galaxies in the fields of the clusters (§4) and found that 
6.5% of IRAC sources had multiple R-band counterparts 
making blending about 1.5 times more common in cluster 
fields. Even though the rate of blends is higher, it should 
not have a significant effect on the LFs. Even in the worst 
case that all 6.5% of IRAC-detected galaxies with multi- 
ple R-band matches are blended (not just coincidentally 
aligned with a faint R-band galaxy in the foreground), 
and those blends are with a galaxy of comparable lumi- 
nosity, the values of M* measured from the LFs would be 
only ~ 0.05 mag brighter. Given the Schechter function 
shape of the LF, it is more probable that most galax- 
ies are blended with a fainter galaxy and therefore 0.05 
mag is likely to be the upper limit of how significantly 
blending affects the LFs. This effect is smaller than the 
statistical errors in the measurement of M* for the LFs 
(§5) and therefore we make no attempt to correct for it, 
but note that our M* values could be systematically high 
by as much as 0.05 mag. 

The large IRAC PSF means that star-galaxy separa- 
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tion using these data is difficult and therefore the clas- 
sification of each matched object is determined from the 
R-band data using the CLASS_STAR parameter from 
SExtractor. This is done using the criteria suggested in 
Fadda et al. (2004). All objects with R < 23.5 with 
CLASS_STAR < 0.9 are considered galaxies. For fainter 
objects with R > 23.5, those with CLASS_STAR < 0.85 
are considered galaxies. Most stars have R - 3.6/um col- 
ors of ^ in the Vega system. The R-band data is 5 
mag deeper than the IRAC data and therefore most stars 
detected by IRAC should be robustly removed using this 
classification. 



2.3. Galaxy Colors 

The most important ingredient in the cluster red se- 
quence algorithm is the measurement of accurate colors. 
Excess noise in the colors causes scatter in the cluster red 
sequence and reduces the probability that a cluster will 
be detected. For images with large differences in seeing, 
PSF shape, and pixel size such as the R-band and 3.6yum, 
measuring accurate colors can be problematic. To this 
end, significant effort was invested in finding the most 
appropriate way to measure colors with this filter com- 
bination. 

Studies of the cluster red sequence using tele- 
scopes/filters with equivalently large angular resolution 
differences (e.g., HST -|- ground based, Holden et al., 
2004; Optical and low-resolution IR, Stanford et al., 
1998) have typically measured colors by degrading the 
highest resolution images using the PSF of the lowest 
resolution images. This is the most accurate way to mea- 
sure colors, and is feasible for a survey of several clusters; 
however, it is time consuming for a survey the size of the 
FLS that has more than a million sources detected in the 
R-band. More importantly, because there are so many 
more galaxies detected in R-band than in 3.6/im, degrad- 
ing those images causes numerous unnecessary blends 
of R-band galaxies resulting in an increased number of 
catastrophic color errors. Degrading the resolution also 
inhibits the potential for detecting distant clusters be- 
cause the signal-to-noise ratio of the faintest R-band ob- 
jects becomes much worse when they are smeared with 
a large PSF. 

The compromise is to use a fixed aperture that pro- 
vides accurate colors, yet is as large as possible for the 
IRAC data (to reduce the need for aperture corrections), 
and yet is as small is possible to reduce the excess sky 
noise in the R-band measurement. It is important to 
use the same diameter apertures for both 3.6/xm and R- 
band so that the colors of bright resolved galaxies are 
measured properly. Galaxies which are small and mostly 
unresolved require an aperture of only 2-3 times the see- 
ing disk to measure a correct color. In principle, colors 
for such galaxies can be measured correctly using a dif- 
ferent sized apertures for both 3.6(Um and R-band (i.e., 
optimized apertures). However, because measuring the 
color correctly for large galaxies that are resolved in both 
filters requires that the aperture must be the same size 
in both filters, we use the same aperture for all galaxies. 

After experimenting with apertures ranging in diame- 
ter between one to ten IRAC pixels (1.22" to 12.2") we 
determined that the three IRAC pixel diameter aperture 
(3.66") was the optimum aperture because it requires a 




Fig. 1. — Color-Magnitude diagram within a 1 Mpc (2.5 arcmin) 
diameter of FLS J171648+5838.6 (cluster 44, Zspec = 0.573), the 
richest cluster in the FLS. Several field galaxies with R - 3.6/im 
colors > 5.5 have been removed for clarity. The solid line is the 
best red sequence model for the cluster (§3.1). The intrinsic scatter 
in the red sequence for this cluster is 0.091 mag for galaxies with 
3.6/im < 17 mag, and 0.151 mag for galaxies with .i.dfim > 17 mag, 
and is comparable to the scatter in other clusters at this redshift. 

relatively small aperture correction at 3.6/im (^ 10%) yet 
is only slightly larger than three R-band seeing disks, re- 
sulting in only a marginal excess sky noise being added to 
the R-band aperture magnitudes. Using this large fixed 
aperture means that the photometry of faintest R-band 
galaxies is not optimized because much of the aperture 
contains sky. As a result, some potential in discovering 
the most distant clusters is sacrificed because the faintest 
red galaxies (i.e., distant red sequence galaxies) may have 
excessively large photometric errors. However, most im- 
portantly, accurate colors are determined for all galaxies, 
and overall the approach provides much better photom- 
etry than degrading the entire survey data. 

As an illustration of the quality of colors achievable 
with this approach we show the color-magnitude diagram 
of FLS J171648-h5838.6, the richest cluster in the survey, 
in Figure 1. The typical intrinsic scatter of early type 
galaxies on the red sequence at the redshift of this cluster 
{zspec = 0.573) is ~ 0.075 ( Stanford et al. 1998; Holden 
et al. 2004). As a comparison we measure the intrinsic 
scatter for FLS J171648-I-5838.6 by subtracting the mean 
photometric error from the total scatter in quadrature. 
This is slightly less rigorous than the Monte-Carlo meth- 
ods used by other authors, but provides a reasonable 
estimate of the scatter. For galaxies with 3.6/im < 17 
mag (3.6/xm > 17 mag) the observed scatter of the red 
sequence is 0.149(0.225) mag, and the mean photometric 
color error is 0.118(0.167) mag, resulting in an intrin- 
sic scatter of 0.091(0.151) mag. We note that without 
knowing the morphologies of the galaxies we are unable 
to properly separate early type galaxies from bluer disk 
galaxies, and therefore this measurement of the scatter is 
almost certainly infiated by Sa or Sb galaxies bluer than 
the red sequence. In particular, these galaxies arc gener- 
ally more prevalent at fainter magnitudes (e.g., EUingson 
et al. 2001). However, even without morphological sep- 
aration, the scatter in the color of red sequence galaxies 
are in fair agreement with scatter in the colors of typical 
red sequences and demonstrates that the 3.66" aperture 
works well for measuring colors. 
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2.4. Keck, WIYN, & SDSS Spectroscopic Data 

A large number of spectroscopic redshifts are available 
for galaxies in the FLS field from several spectroscopic 
campaigns. A sample of 642 redshifts were obtained us- 
ing the HYDRA spectrograph on the Wisconsin-Illinois- 
Yale-NOAO (WIYN) 3.6m telescope as part of a pro- 
gram to foUowup radio sources in the FLS (Marleau et 
al. 2007). A set of 1373 redshifts in the FLS field were ob- 
tained for galaxies selected by their red R - Kg colors us- 
ing the DEIMOS spectrograph on the 10m KECK II tele- 
scope by Choi et al. (2006). Lastly, 1296 redshifts were 
obtained using the Hectospec Fiber Spectrograph on the 
6.5m MMT by Papovich et al. (2006). The primary tar- 
get of that survey were galaxies that are detected in the 
FLS MIPS 24/im imaging and have R < 21.5 mag. In 
addition to redshifts from these projects, 1192 redshifts 
in the FLS field are also available from the Sloan Digital 
Sky Survey (SDSS) DR5 database (Adelman-McCarthy 
et al. 2007). In total there are 4503 redshifts at various 
positions available in the FLS. Of these, 26 are likely to 
be cluster red sequence galaxies (see §3.7). 

2.5. Palomar Spectroscopy 

In addition to the spectroscopic catalogues available, 
we also obtained our own longslit spectroscopy for bright 
red sequence galaxies in three clusters with 0.4 < Zphot < 
0.6 in the FLS using the Double-Spectrograph (Double- 
spec) on the 200-inch Hale telescope at Palomar Moun- 
tain (P200). We also obtained multi-object spectroscopy 
using the COSMIC Spectrograph on the P200 for an ad- 
ditional three clusters with Zphot < 0.3. These six clusters 
were chosen for foUowup because they were amongst the 
richest clusters in our preliminary cluster catalogues. 

2.5. L Double-Spectrograph Data 

Spectroscopy of bright red sequence galaxies in clusters 
FLS J17124H-5855.9, FLS J172122-h5922.7, and FLS 
J171648-h5838.6 (clusters 16, 38, and 44 listed in Table 1) 
was performed on 17, 18, and 19 August 2004 with Dou- 
blespec on the P200. The observations were made with 
the "Red" camera using the 316 1/mm grating blazed at 
7150A and a 0.5" wide slit, giving a spectral resolution 
of 2.6A (^ 150 km s~^). The Doublespec longslit is ^ 
1.5 arcmin long and the angle of the slit on the sky can 
be rotated. In all 3 clusters we centered the slit on the 
brightest cluster galaxy (BCG) and then chose a rotation 
angle so that we could get at least 2 other bright objects 
(preferentially red sequence galaxies) on the slit. 

For FLS J172122-h5922.7 and FLS J171648-I-5838.6 
we obtained spectra of 3 objects in the field, and in 
FLS J17124H-5855.9 we managed 4. We obtained three 
20 minute exposures for FLS J172122-)-5922.7 and FLS 
J171648-I-5838.6, which have photometric redshifts of 
0.57 and 0.55, respectively, and one 20 minute exposure 
for FLS J171241-t-5855.9, which has a photometric red- 
shift of 0.39. We also observed a spectroscopic standard, 
calibration lamps, dome flats and twilight flats at the 
beginning of each night. Data reduction and wavelength 
calibration were performed using the standard IRAF 
techniques. After 1-d spectra were extracted, 7 of the 
10 objects had a signal-to-noise ratio suitable for cross- 
correlation. One of the spectra in FLS J171241-I-5855.9 
has a strong emission line at 7056A and no possible iden- 



tification that puts it near the photo-z of the cluster. 
This object was therefore considered a field interloper. 
The remaining 6 spectra (two per cluster) showed signif- 
icant absorption features typical of early type galaxies 
and redshifts were obtained by cross correlating them 
with an elliptical galaxy spectrum. The redshifts of the 
galaxies within each cluster were similar (Az < 0.01) and 
are in excellent agreement with the cluster photometric 
redshift. These spectroscopic redshifts are listed in the 
cluster catalogue (Table 1). 

2.5.2. COSMIC Data 

Multi-object spectroscopy of both red sequence galax- 
ies and MIPS 24/im-detected galaxies in the fields of clus- 
ters FLS J171059-h5934.2, FLS J171639-H5915.2, FLS 
J171505-h5859.6, and FLS J172449-H5921.3 (clusters 1, 
2, 8, and 10 listed in Table 1) were performed on 26, 
27, 28, 29 May 2006, and 15, 16, 17 June 2007 using 
the COSMIC Spectrograph on the 200" Hale Telescope 
at Palomar Mountain. These observations were made 
with the 300 1/mm grating blazed at 5500A with 1" 
wide slits giving a spectral resolution of 8A (~ 450 km 
s~^). These data are part of a larger campaign to study 
cluster 24/im sources and full details of the data reduc- 
tion, calibration and cross-correlation will be presented 
in a future paper (Muzzin et al. 2008, in preparation). 
We obtained 17, 16, 12 and 20 good-quality spectra in 
the fields of FLS J171059-t-5934.2, FLS J171505-I-5859.6, 
and FLS J172449+5921. 3 respectively, and redshifts were 
determined using cross-correlation. Including the data 
from the other spectroscopic campaigns, the field of FLS 
J171059-H5934.2 has 10 galaxies with z = 0.126, the field 
of FLS J171639+5915.2 has 7 galaxies with z = 0.129, 
the field of FLS J171505+5859.6 has 9 galaxies with z 
= 0.252, and the field of FLS J172449-f 5921.3 has 12 
galaxies with z = 0.253. These redshifts are included in 
the cluster catalogue (Table 1). 

2.6. Keck/DEIMOS Spectroscopy of FLS 
J172126+5856.6 

Spectroscopy was obtained of the candidate cluster 
FLS J172126-I-5856.6 (cluster 93 in Table 1) with the 
Deep Imaging Multi-Object Spectrograph (DEIMOS, 
Faber et al. 2003) on the 10 m Keck II telescope. On the 
night of 1 September 2005, we obtained three 1800s ex- 
posures on the same mask in non-photometric conditions 
with -1.3" seeing. The 600ZD grating (Abiazo = 7500 
A; AAfwhm — 3.7 A) and a GG455 order-blocking fil- 
ter were used. The DEIMOS data were processed using 
a slightly modified version of the pipeline developed by 
the DEEP2 team at UC-Berkeley^^. Relative flux cali- 
bration was achieved from observations of standard stars 
from Massey & Gronwall (1990). 

Slits were preferentially placed on candidate red se- 
quence galaxies, allowing a total of 10 slits on likely clus- 
ter members. Of the 10 candidate red sequence galaxies, 
five had sufficient S/N for determining redshifts, and four 
had redshifts Az < 0.01 from each other, with the z = 
1.045. These redshifts are included in the cluster cata- 
logue (Table 1). 

|http : //astro ■ berkeley . edu~cooper/deep/spec2d/ 1 
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3. CLUSTER FINDING ALGORITHM 

The cluster finding algorithm employed in this study is 

essentially the CRS algorithm of Gladders & Yee (2000, 
2005) with some minor modifications. Here we outline 
only the major steps, and refer to those papers for a 
more detailed explanation of the procedures. 

The CRS algorithm is motivated by the observation 
that early type galaxies dominate the bright end of the 
cluster LF and that these galaxies always follow a tight 
red sequence in the color-magnitude plane. At increas- 
ing redshift the observed red sequence color becomes red- 
der^^ and because this change in color follows closely the 
predictions from a passively evolving stellar population, 
the color can be used as a robust photometric redshift 
estimate for a cluster. In order to apply the CRS algo- 
rithm, slices are made in the color-magnitude plane of a 
survey. Galaxies are then assigned weights based on the 
probability that they belong to a particular slice. This 
probability is determined by the color and the photomet- 
ric error in the color. Once color weights for each galaxy 
in each slice have been assigned, each galaxy is also as- 
signed a magnitude weight. Magnitude weighting is done 
because bright red sequence galaxies are more likely to 
be members of clusters than faint ones. 

Once each galaxy is assigned a color and magnitude 
weight for each slice, the positions of each galaxy are plot- 
ted for each slice with their respective weights. The re- 
sulting "probability map" for each slice is then smoothed 
and peaks in these maps represent likely cluster candi- 
dates. In the following subsections we discuss in more 
detail the steps in our version of the algorithm. 

3.1. Red Sequence Models 

The first step in finding clusters with the CRS is to 
model the color, slope, and intercept of the cluster red 
sequence as a function of redshift. This was done by 
making simulated single-burst galaxies using all available 
metallicities from the Bruzual & Chariot (2003) spectral 
synthesis code. The models are constructed with 50% 
of the stars forming in a single-burst at t = 0, and the 
remainder forming with an exponentially declining star 
formation rate of r = 0.1 Gyr. Using a range of metal- 
licities causes the color of each galaxy to be slightly dif- 
ferent at z = 0, with the most metal-rich galaxies being 
the reddest. The absolute magnitude of each galaxy with 
a different metallicity is normalized using the U-V, V-I, 
and J-K red sequences of Coma (Bower et al. 1992) as- 
suming that a metallicity gradient with magnitude is the 
primary source of the slope of the red sequence. Nor- 
malizing the absolute magnitude of each galaxy this way 
allows us to reproduce models with the correct red se- 
quence color and slope with redshift. 

There is increasing evidence that the slope of the red 
sequence is not only caused by a metallicity scqiience, but 
is also the prochict of an age sequence, with the less Imni- 
nous galaxies being both more metal poor and younger 

^'^ The observed-frame color of the red sequence becomes redder 
with increasing redshift because of band shifting. Tlic rest-frame 
color change due to passive evolution actually makes galaxies bluer 
at higher redshift, but is a small effect for a single-burst popula- 
tion formed at high redshift. Because the change in observed-frame 
color is dominated by the k-corrcction from an old stellar popula- 
tion, it increases monotonically with redshift and provides a good 
estimate of the cluster redshift. 



(e.g., Nelan et al. 2005; Gallazzi et al. 2006). Examina- 
tion of spcctroscopically confirmed clusters in the FLS 
shows that the pure metallicity sequence used in our 
models reproduces the red sequence slope and color quite 
well, and because we are only interested in a fiducial red 
sequence model for detecting clusters and determining 
photometric redshifts, no further tuning of the ages of 
galaxies along the sequence is done. 

Once the absolute magnitude of each model galaxy is 
normalized using the Coma red sequences, linear fits to 
the R - 3.6^m vs. 3.6/xm color-magnitude relations of 
the model galaxies between 0.1 < z < 1.6 are made. A 
high density of redshift models are fit so that there is 
significant overlap in color space (185 slices between 0.1 
< z < 1.6). This assures that no clusters are missed 
because they have colors between the finite number of 
models and it also allows for increased precision in the 
photometric redshifts. 

We computed two sets of single-burst models, one with 
a formation redshift {zf) = 2.8, and another with Zf = 
5.0. These two sets of models produce nearly identical 
observed red sequences at z < 1.1, but begin to diverge at 
higher redshifts. There is evidence from previous studies 
of the fundamental plane (e.g., van Dokkum et al. 1998; 
van Dokkum & Stanford 2003, and many others), evo- 
lution of the color-magnitude diagram (Stanford et al. 
1998; Holden et al. 2004), and K-band luminosity func- 
tion (De Propris et al. 1999; Lin et al. 2006; Muzzin et 
al. 2007a) that a ~ 3 model is appropriate for cluster 
early types; however, the uncertainties in these studies 
are fairly large. There is also evidence that many of the 
most massive field early types formed the majority of 
their stars at z > 5 (McCarthy et al. 2004; Glazebrook 
et al. 2004), so the possibility remains that a Zf = 5.0 is 
more appropriate. Regardless, the majority of the sys- 
tems we have discovered are at z < 1.1, and therefore the 
Zf uncertainty does not affect the photometric redshifts 
of these systems. For systems at z> 1.1, the redshift can 
be considered an upper limit. For example, the photo-z 
for a cluster at z = 1.3 in the Zf = 2.8 model would be z 
~ 1.2 in the Zf = 5.0 model. 

To illustrate the depth of the survey, and the loca- 
tion of the red sequence models. Figure 2 shows the R - 
3.6/im vs. 3.6/im color-magnitude diagram for all galax- 
ies in the FLS with some of the Zf = 2.8 red sequence 
models overlaid. The density of galaxies with M ~ M* 
begins to drop off significantly for the z > 1.2 red se- 
quence models because of the depth of the R-band data 
(M*(3.6/J,m) ^ 17.0 mag at z = 1.2) and therefore we 
consider z ~ 1.2 the upper limit at which we can reliably 
detect clusters. Remarkably, the red sequence models are 
even well separated in color space at z < 0.5 where the 
R - 3.6/xm filters do not span the 4000A break. This is 
caused by the large wavelength separation between the 
bands, and the wide 3.6/xm filter which has a strongly 
redshift-dependent negative k-correction. Although the 
k-correction in R-band evolves slowly with redshift out 
to z ~ 0.5, the k-correction for 3.6/um is significant and 
therefore the R - 3.6/im color is still a sensitive redshift 
indicator. 

3.2. Color Weights 

Once red sequence models have been made, weights 
based on the probability that a galaxy belongs within 
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Fig. 2. — Observed color-magnitude diagram for all galaxies in the FLS. The solid lines are fiducial red sequence models at different 
redshifts generated using the Bruzual & Charot code. The redshift of each model is labelled in the figure. The bulk of the shift in color 

with redshift of the models is due to bandpass shifting or "k-correction" , not because of evolution in the rest-frame colors of the galaxies. 



a color slice are computed. The typical la scatter in 
the local cluster color-magnitude relation is ~ 0.075 mag 
(e.g., Lopez-Cruz et al. 2004; Bower et al. 1992). The 
scatter has been measured in clusters to z ~ 1 where it 
remains remarkable consistent (e.g., Stanford ct al. 1998; 
Gladders et al. 1998; Blakeslee et al. 2003). Thereafter, 
it may become somewhat more scattered (Holden et al. 
2004). Assuming that this relation holds to ^ ~ 1.3, color 
weights (with values ranging from to 1) are assigned by 
computing the overlapping area of a galaxy's color with 
the red sequence assuming a red sequence intrinsic dis- 
persion of 0.075 mag and assuming the galaxy's color 
is represented by a Gaussian centered on the measured 
color with a Icr dispersion equal to the color error (see 
e.g. GYOO, Figure 3 for an example). Using this method, 
the weight of a bright galaxy lying directly on the red se- 
quence with a color error significantly narrower than the 
width of the red sequence is 1.0. The same galaxy, with 
a color error equal to the dispersion in the red sequence, 
has a weight of 0.67. Color weights are computed for all 
galaxies in all 185 color slices. 

3.3. Magnitude Weights 

In addition to the color weights, galaxies are also 
weighted based on their magnitude rc;lative to a fiducial 
M* value. Cluster early types are usually the bright- 
est, reddest galaxies at a given redshift and therefore, 
the brightest galaxies within a color slice are more likely 
to be cluster galaxies and should be given extra weight. 
The distribution of magnitude weights was defined as 



P(M) by GYOO (see their §4.3). We compute the P(M) 
using the data themselves, as suggested by those au- 
thors, and when doing so we consider objects within the 
one-percentile highest density regime as "cluster" galax- 
ies. This is a slightly more strict cut than the ten- 
percentile cut used by GYOO; however, the fact that 
IR-selected galaxies are more strongly clustered than 
optically-selected galaxies justifies using a more stringent 
cut. 

3.4. Probability Maps 

Once the magnitude and color weights for all galaxies 
in each of the individual color slices have been computed, 
a probablility map of each slice is created. The map is 
a spatial galaxy density map of the survey within each 
redshift slice. The map is made using pixels which are 
125 kpc in physical size at the redshift of each slice. The 
probablility flux from each pixel is determined by placing 
each galaxy on the pixel that corresponds to its location 
in the survey, weighted by the product of its color and 
magnitude weights. Once each slice is constructed this 
way, it is smoothed with the exponential kernel suggested 
in GYOO (their equation 3). 

3.5. Noise Maps 

The noise properties of the probability maps of differ- 
ent color /redshifts slices are usually different. In partic- 
ular, the maps of the highest redshift slices tend to have 
large noise peaks because the survey is only as deep as ~ 
M* in those slices. The lower redshift probability maps 
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have a smoother background because there arc numer- 
ous M > M* galaxies which are more evenly distributed 
spatially and have a low probability of belonging to any 
slice because they have a large color error. The higher 
redshift maps are shallower, thereby lacking the M > 
M* galaxies which provide this smooth background. If 
peak finding is run on all probability maps using similar 
detection parameters, it produces significantly different 
numbers of detections in different slices. In particular, 
almost any noise in the highest redshift maps results in 
the detection of a "cluster" . 

To circumvent this problem the parameters of the peak 
finding for each map can be tuned individually in order to 
produce a reasonable number of detections in each slice; 
however, the resulting cluster catalogue is clearly biased 
by what is considered a "real" detection in a given map. 
It is preferrablc to have a cluster catalogue which is as 
homogenously-selected as possible and based on a quanti- 
tative selection. Therefore "noise" maps are constructed 
and are added to each probability map to homogenize 
their noise properties. 

The noise maps arc constructed by adding fake red se- 
quence galaxies to each pixel of the probability maps. 
Adding a constant background of fake galaxies does not 
change the noise properties of a map because it is the 
variance in the number of background galaxies that de- 
termines the noise. Wc experimented with a variety of 
variances to add, but settled that the variance from the 
photometric color errors of six M* red sequence galaxies 
per pixel provided the best results. This level of noise 
removes the spurious detections in the highest redshift 
slices, but docs not add so much noise as to wash-out 
the majority of the poorer clusters in the lower redshift 
slices. 

The noise in each pixel is calculated by first deter- 
mining the average color error of an M* red sequence 
galax;y using the survey data. Once the average color 
error per slice is tabulated, the weights of six M* red se- 
quence galaxies are Monte Carlo simulated for each pixel 
of a noise map assuming that the colors are normally 
distributed around the red sequence with a dispersion 
equal to the mean color error. These simulated weights 
are then assigned to each pixel of the noise map and 
each noise map is added to the appropriate cluster prob- 
ability map. This approach thereby implictly defines a 
"cluster" as an overdensity detectable above the Poisson 
noise from six M* background red sequence galaxies at 
any redshift. The noisc+clusters maps have similar noise 
properties for every slice and peak finding can be run us- 
ing identical parameters for all maps. 

We note that in our simple empirical method for ho- 
mogenizing the noise in the probability maps the added 
noise is Poissionian, not clustered like the underlying 
background galaxy distribution. Despite this, the noise 
maps technique works extremely well, effectively smooth- 
ing out spurious noise spikes in the highest redshift prob- 
ability maps. In principle, a more sophisticated method 
which includes the clustering properties of background 
galaxies could be implemented; however, for our pur- 
poses such an approach is unecessary. Only the detection 
probability of poorest clusters near the significance limit 
are affected by different choices in noise maps. Galaxies 
from the poorest clusters do not contribute significantly 
to the LPs, which are dominated by counts from more 



massive systems, and therefore we do not consider this 
issue further. 

3.6. Cluster Detection 

Once the combined noise-probability maps have been 
made, peaks are detected in each map using SExtractor. 
The peak-finding is done differently from GYOO in that 
the individual 2d slices are searched instead of merging 
the slices into a 3d datacube and searching for 3d peaks. 
It is unclear how these two methods compare; however, 
they are likely to be similar and searching the slices in- 
dividually permits easy visual inspection of the sources 
on each map which allows us to check any problems that 
have occured with peak finding or in the generation of 
the map. Pixels 5a above the background are fiagged 
and 25 connected pixels are required to make a detec- 
tion. 

The slices are close in color space and therefore clus- 
ters (particularly rich ones) are detected in more than 
one color slice. The same cluster is identified in mul- 
tiple color slices by merging the slice catalogues using 
a matching radius of 8 pixels (1 Mpc). Clusters found 
across as many as 20 color slices arc connected as being 
the same object. The color slices are not linear in red- 
shift, but 20 slices correspond to Az ^ ±0.06. These 
combined spatial and color limits for connecting clusters 
imply that clusters with separations > 1 Mpc in trans- 
verse distance and > 0.06 in redshift space can be re- 
solved into distinct sy stems This level of sensitivity is 
similar to that found by Gladders & Yec (2005) using R - 
z' colors to select clusters. They also demonstrated that 
subclumps at redshift spacings much less than this are 
likely to be associated subclumps or infalling structures 
related to the main body of the cluster. 

3.7. Photometric Redshifts 

Each cluster is assigned the photometric redshift of 
the color slice in which it is most strongly detected. The 
strength of the detection is determined by using SEx- 
tractor to perform aperture photometry of each cluster 
on each probability map. This provides a "probability 
flux" , and the cluster is assigned to the slice in which it 
has the largest probability flux. 

The large number of spectroscopic redshifts available 
for the FLS can be used to verify the accuracy of the red 
sequence photometric redshifts. Examining the spectro- 
scopic catalogue for galaxies within a 1 Mpc circle around 
each cluster shows that there are numerous galaxies with 
spectroscopic redshifts in the field of many of the clusters. 
The spectroscopic targets were chosen with a variety of 
selection criteria (none of which preferentially select early 
type galaxies) and therefore the majority of galaxies with 
redshifts are foreground or background galaxies. We use 
only the spectroscopic redshifts for galaxies which have 
a combined magnitude and color weight of > 0.2 in or- 
der to preferentially select likely cluster members. This 
cut in weight is used because it corresponds to the typ- 
ical combined magnitude and color weight of M < M* 

The color slices are closer together at 2: > 1 and only systems 
with Az > 0.12 can be resolved at this redshift. We note that 
although the overall level of projections is likely to be low, because 
of the bunching up of the color slices, the highest redshift clusters 
will be the most suseptible to projection effects. 
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Fig. 3. — Photometric vs. spectroscopic redshift for clusters 
in the FLS field. The color of the circle corresponds to the tele- 
scope/project where the spectroscopic redshifts were obtained (see 
§2). Large circles denote clusters with more than one spectroscopic 
redshift, small circles denote clusters with only one spectroscopic 
redshift. Excluding the one large outlier at Zspec ~ 0.9, the rms 
scatter is Az = 0.04. 

red sequence galaxies. Once the cut is made there are 23 
clusters which have at least one spectroscopic redshift for 
a likely cluster red sequence galaxy. Remarkably, there 
are 26 galaxies which meet this criteria and 24 of these 
have a spectroscopic redshift < 0.1 from the photometric 
redshift of the cluster. This illustrates the effectiveness 
of the red sequence color at estimating photometric red- 
shifts provided that the galaxy has a high-probability of 
being a cluster early type. 

In Figure 3 we plot spectroscopic vs. photometric red- 
shift for these 23 clusters plus the additional 6 for which 
we obtained our own spectroscopic redshifts (§ 2.5). The 
straight line marks a one-to-one correlation. Large points 
represent clusters with more than one galaxy with a red- 
shift consistent with the being in the cluster. Small 
points represent clusters with a single spectroscopic red- 
shift. Excluding the large single outlier with z^pec ~ 0.9 
(which is likely to be a bluer galaxy at high-redshift based 
on its spectrum and IRAC colors, see §6.3) the rms scat- 
ter in the cluster spectroscopic vs. photometric redshift 
is Az = 0.04, demonstrating that the photometric red- 
shifts from the red sequence algorithm work extremely 
well. 

The accuracy of the photometric redshifts from the 
FLS sample is comparable to the accuracy of the RCS 
surveys (Yee et al. 2007; Gladders & Yee 2005) which 
use R - z' color selection, even though the R - 3.6/im col- 
ors have larger photometric errors than the R - z' colors. 
It is likely this is because the model red sequence colors 
change much more rapidly with redshift in R - 3.6/im 
than in R - z' (R - 3.6/im spans 2 magnitudes between 
z = 0.5 and z — 1.0, whereas it spans only 1 magnitude 
in R - z'). The larger change in the R - 3.6/im colors 
with redshift means that photometric measurement er- 
rors should correspond to smaller errors in photometric 
redshift. 

3.8. Bgc Richness Parameter 

The final step in the selection of the cluster sample is to 
cut low-richness detections from the catalogue. The false 
postive rate is higher for low richness systems (i.e., galaxy 
groups) and we prefer restrict our analysis of the cluster 
LFs to a high confidence sample of massive clusters. The 



cluster richnesses are measured quantitatively using the 
Bgc richness parameter (Longair & Seldner 1979; for a de- 
tailed look at the application of Bgc to measuring cluster 
richnesses see Yee & Lopez-Cruz 1999). Bgc is the am- 
plitude of the 3-dimensional correlation function between 
cluster galaxies and the cluster center. Bgc is measured 
within a fixed aperture (typically 500 kpc radius) and is 
well-correlated with cluster physical parameters such as 
velocity disperion (cr). X-ray temperature (T^^), and the 
radius as which the mean density of the cluster exceeds 
the critical density by a factor of 200 (R200) (e.g., Yee & 
Lopez-Cruz 1999; Yee & EUingson 2003; Gilbank et al. 
2004; Muzzin et al 2007b). 

Gladders & Yee (2005) introduced a new form of the 
Bgc parameter, counting the overdensity of red sequence 
galaxies within a fixed aperture, rather than all galax- 
ies and defined this new parameter as Bgc,_R. This form 
of richness suffers less from cosmic variance in the back- 
ground because red sequence galaxies provide better con- 
trast with the field, and therefore it is a more robust es- 
timate of the cluster richness. We use Bgc^R rather than 
Bgc for determing the richnesses of the FLS clusters. The 
net number of 3.6/xm red sequence galaxies with M < M* 
+ 1.0 (where M* is determined from the data itself, see 
§5.1) are counted within a fixed aperture of 500 kpc ra- 
dius. The model red sequences from §3.1 are used and 
galaxies within ± 0.3 in color are considered to belong 
to the red sequence. 

Systems with Bgc,_R < 200 are removed from the clus- 
ter catalogue. The Bgc-M2oo relation of Yee & EUingson 
(2003) implies that this corresponds to removing groups 
with M200 < 6.6 X 10^'^ Mq, where M200 is defined as 
the mass contained with R2oo- Groups with masses be- 
low this typically have only ~ 5-10 bright galaxies (e.g., 
Balogh et al. 2007), making them difficult to select ro- 
bustly with the CRS algorithm. The systems that are 
removed by the richness cut are typically tight compact 
groups of 3-4 extremely bright galaxies that are the same 
color. Although they are not rich, they have a strong 
probability of being detected by the CRS algorithm be- 
cause of their luminosity and compactness. It is likely 
that the majority of these systems are bona fide low- 
richness galaxy groups; however, we have no way of ver- 
ifying the false-positive rate for these systems. 

Before these low-richness systems are cut from the 
catalogue there are 134 cluster candidates between 0.1 
< z < 1.4 in the FLS field. Removing systems with 
Bgcfl < 200 leaves a total of 99 candidate clusters in the 
sample. 

3.9. Cluster Centroids 

Defining a centroid for clusters can be a challenging 
task, yet is extremely important because properties de- 
termined within some aperture around the cluster (such 
as richness, or LF) can vary strongly with the choice of 
centroid. In many cluster studies the location of the BCG 
is used as the center of the cluster. This is a reasonble 
definition as frequently the BCG lies at the center of the 
dark matter halo and X-ray emission; however, there are 
also many examples where it does not. Furthermore, not 
all clusters have an obvious BCG, particularly at higher 
redshift. 

Given these issues, two centroids are computed for the 
FLS clusters, one based on the location of the peak of 
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the red sequence probability flux in the probabiUty maps, 
and the other based on the location of the BCG within 
500 kpc of this centroid. In order to avoid bright fore- 
ground galaxies the brightest galaxy in the field with a 
red sequence weight > 0.4 is designated as the BCG. 
Eye examination of the clusters shows that this criteria 
is effective at choosing what appears visually to be the 
correct galaxy; however, because it chooses only a single 
galaxy this technique is still potentially suseptible to red 
low-redshift field interlopers. 

When computing the cluster LPs, only one of the cen- 
troids can be used. We define an "optimmn" centroid for 
each cluster using the ^gc,R parameter. Bgc.R is com- 
puted at both centroids and the optimum centroid is the 
centroid which produces the max;imum value of Bgc,ij. 
This approach is simplistic, but because Hgc is the corre- 
lation amplitude between the cluster center and galaxies, 
the centroid which produces the largest value should be 
the best centroid of the galaxy population. 

4. PROPERTIES OF THE CLUSTER CATALOGUE 

The final cluster catalogue of 99 clusters and groups is 
presented in Table 1. Where spectroscopic redshifts are 
available for high-probability cluster members they are 
listed in column 3, with the number of redshifts in paren- 
thesis. For each cluster we also compute an estimate of 
R200 and M200. The M200 values are estimated using the 
correlation between Bgc and M200 measured by Muzzin 
et al. (2007b) for 15 X-ray selected clusters at 2; ^ 0.3 in 
the K-band. The K-band and 3.6/tim bandpasses sample 
similar parts of a galaxy's spectrum at 0.1 < z < 1.5 
and therefore it is reasonable to assume that Bgc values 
measured in both these bands will be comparable. The 
best-fit relation between M200 and Bgc is 

Log{M20Q) = (1.62±0.24)Lo5(Bgc) + (9.86 ±0.77). (1) 

Muzzin et al. (2007b) did not measure the correlation 
between Bgc and R200 in the K-band, although Yee & 
EUingson (2003) showed a tight correlation for the same 
clusters using r-band selected Bgc. Using the Muzzin et 
al. (2007b) K-band data we fit the correlation between 
these parameters for those clusters and find that the best 
fit relation is 

LogiR2oo) = (0.53 ± 0.09)LogiBgc) - (1.42 ± 0.29). (2) 

The rms scatter in the M200 - Bgc relation is 35% and 
for the R200 - Bgc relation it is 12%. These scatters are 
similar to that measured between M200 and K-band se- 
lected richness (parameterized by N200) at z ~ by Lin 
et al. (2004). The values of M200 and R200 derived from 
these equations are listed in columns 9 and 10 of Table 
1, respectively. 

We caution that these equations have only been cal- 
ibrated using rich clusters, and that extrapolating to 
lower richness clusters such as those in the FLS may not 
be appropriate. The lowest richness cluster in the Muzzin 
et al. (2007b) sample has a richness of Log(Bgc) = 2.8, 
yet the majority of clusters in the FLS (70/99) have lower 
richnesses than this. There is evidence from both obser- 
vations (e.g., Lin et al. 2004) and numerical simulations 
(e.g., Kravtsov et al. 2004) that the same power-law cor- 
relation between cluster galaxy counts (which are closely 
related to Bgc) and M200 extend to richnesses well lower 
than our Bgc,ij > 200 cut, and therefore it probably not 



too unreasonable to extrapolate equations (1) and (2) to 
lower richnesses. 

Using an indirect method to estimate M200 and R200 
means that reliable errors in R200 and M200 can not be 
computed for individual clusters; however, the rms scat- 
ters in the correlations are at least indicative of the av- 
erage uncertainty in the measurement of the parameters 
for the sample. Therefore, we suggest that the average 
error in the M200 and R200 values listed in Table 1 are ± 
35% and 12% respectively, but that the error in a par- 
ticular cluster can be several times larger or smaller. 

In Figure 4 we plot a histogram of the number of clus- 
ters as a function of redshift in the FLS. The solid his- 
togram shows the distribution of all clusters and the dot- 
dashed histogram shows the distribution of clusters with 
M200 > 3 X 10^'^ Mq (BgcR > 700). Similar to the predic- 
tions of numerical simulations (e.g., Haiman et al. 2003) 
the number of clusters peaks at z ~ 0.6. Qualitatively, 
the distribution of clusters is also similar to that found 
by Gladders & Yee (2005) in comparable size patches; 
however, the cosmic variance in the number of clusters 
in 4 deg^ patches is too large to make a meaningful 
comparison between the selection of clusters in the R - 
z' bandpasses versus the R - 3.6y^m bandpasses. 

We plot the locations of the clusters superposed on the 
3.6/j.m image of the FLS field in Figure 5 as open circles. 
Large and small circles represent clusters with M200 > 3 
X lO^"* Mq and M200 < 3 x 10^'' Mq, respectively, and 
clusters with photometric redshifts 0.1 < 2; < 0.4, 0.4 
< z < 0.8, z > 0.8 are plotted as blue, green, and red 
circles, respectively. The clusters themselves are clearly 
clustered; demonstrating the need for wide-field surveys 
when searching for representative samples of galaxy clus- 
ters. 

We show a few examples of some of the richest cluster 
candidates in Figures 6-11. The top left panel for each 

Figure is the R-band image; the top right is the 3.6/xm 
image, and the bottom left panel is the 8.0/xm image. 
All images are 1 Mpc across at the; cluster redshift. The 
bottom right panel of each figure shows a histogram of 
the color distribution of galaxies with M < M* within a 
1 Mpc diameter aperture. The color of the red sequence 
for the photometric redshift is marked with an arrow. 
The dashed histogram is the mean color background in 
that aperture measured from the entire survey. The er- 
ror bars on the dashed histogram are computed as the la 
variance in each bin from 200 randomly selected 1 Mpc 
apertures within the survey. Galaxies are clustered, and 
therefore assuming the variance is Gaussian-distributed 
is probably an overestimate of the true variance (because 
there will be large wings in the distribution due to clus- 
tering); however, it provides a first-order demonstration 
of the overdensity of the cluster relative to the field. 

Overall, the cluster catalogue is qualitatively similar 
in both redshift, and richness distributions to catalogues 
selected with the same technique in diff'erent bandpasses 
(e.g. Gladders & Yee 2005, Gilbank et al. 2004), demon- 
strating that clusters can be reliably selected with the 
CRS method on IRAC data despite the limited spatial 
resolution of the instnmient. 

5. CLUSTER LUMINOSITY FUNCTIONS 

In this section we measure the IRAC luminosity func- 
tions of the FLS cluster sample and use these to study 
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Fig. 4. — Redshift distribution of clusters in the FLS. The solid 
histogram is for all clusters and the dot-dashed histogram is for 

clusters with M200 > 3 x lO^"' Mq. 

the evolution of stellar mass assembly and dusty star for- 
mation/AGN activity in clusters. 

5.1. The 3.6^m and 4 -5^171 Luminosity Functions 

The luminosity of galaxies at 3.6/im and 4.5/im over 
the redshift range 0.1 < z < 1.5 is dominated by emis- 
sion from low mass stars and is fairly insensitive to on- 
going star formation or dust. Consequently, the 3.6/im 
and 4.5/xm cluster LPs provide an estimate of the stellar 
mass function of cluster galaxies, and their redshift evo- 
lution can constrain the mass assembly history of cluster 
galaxies. One concern with using these LFs as a proxy 
for the stellar mass function is that at ^ < 0.5 the 3.3/im 
PAH feature found in strongly star forming galaxies can 
contaminate the stellar emission observed at 3.6/im and 
4.5/im; however, it is likely that such contamination will 
be small for cluster galaxies in this redshift range. In a 
study of Luminous Infrared Galaxies (LIRGs, L/ij > 10^^ 
L0) with estimated star formation rates of ~ 100 Mq 
yr~^, Magnelli et al. (2008) found that the excess emis- 
sion in the IRAC bands due to the 3.3/tm PAH feature 
was only ^ 30%. Given that such luminous LIRGs are 
fairly rare at z < 0.5 (e.g., Perez-Gonzalez et al. 2005), 
and the increase in flux is small, even for strongly star 
forming galaxies, contamination of the 3.6/im and 4.5/im 
bandpasses by 3.3/im PAH emission should be negligible. 

Another concern is that in the worse cases there can 
be variations in the stellar mass-to- light ratio (M*/L) of 
galaxies in similar bandpasses as large as a factor of 5- 
7 (such as in the K-band, e.g., Brinchmann 1999; Bell 
& de Jong, 2001; Beh et al. 2003). These variations 
are smaller for evolved populations such as those found 
in clusters and in general the luminosity of a galaxy at 
3.6/tm and 4.5/im is still a reasonable proxy for its stellar 
mass. 

Exhaustive studies of both the K-band (e.g., De Pro- 
pris et al. 1999, Lin et al. 2006, Muzzin et al. 2007a) 
and 3.6/tm &; 4.5/tm (Andreon 2006, De Propris et al. 
2007) LFs of cluster galaxies have shown that the evolu- 
tion of M* in these bands is consistent with a passively 
evolving stellar population formed at high-redshift (zj > 
1.5), suggesting that the majority of the stellar mass in 
bright cluster galaxies is already assembled into massive 
galaxies by at least 0^1. Here we compute the LFs 
in the 3.6/tm and 4.5/tm bands for the FLS clusters to 
confirm that the FLS cluster sample provides similar re- 
sults, and to demonstrate that these LFs can be used to 



estimate the stellar contribution to the MIR cluster LFs 
(§5.2). 

The LFs are measured by stacking clusters in redshift 
bins of Az = 0.1 starting from z = 0.1. For each clus- 
ter, the number of galaxies within R200 in 0.25 mag bins 
is tabulated and the expected number of background 
galaxies within R200 is subtracted from these counts. 
The background counts are determined from the entire 
3.8 deg^ survey area and are well constrained. Each 
background-subtracted cluster LF is then "redshifted" 
to the mean redshift of the bin using a differential dis- 
tance modulus and a differential k-corrcction determined 
from the single-burst model (§3.1). At 3.6/im and 4.5/im 
the k-corrections for galaxies are almost independent of 
spectral- type (e.g., Huang et al. 2007a) and therefore 
using only the single-burst k-correction rather than a k- 
correction based on spectral-type does not affect the LFs. 
Furthermore, the differential k-corrections and distance 
moduli are small (typically < 0.1 mag) and do not affect 
the LFs in a significant way. 

The final stacked LFs are constructed by summing the 
individual LFs within each bin. The errors for each mag- 
nitude bin of the final LF are computed by adding the 
Poisson error of the total cluster counts to the Poisson 
error of the total background counts in quadrature. 

In Figure 12 and Figure 13 we plot the 3.6/im and 
4.5/im cluster LFs respectively. The 3.6/im LFs are fit 
to a Schechter (1976) function of the form 

ct>{M) = (0.41nlO)<^.(10°-4(^*-^))^+"exp(-10°-^(^*-^)), 

(3) 

where; a is the faint-end slope; (f>^, the normalization; and 
M* is the "characteristic" magnitude, which indicates the 
transition between the power-law behavior of the faint- 
end and the exponential behavior of the bright end. The 
functions are fit using the Levenberg-Marquardt algo- 
rithm for least-squares (Press et al. 1992) and errors are 
estimated from the fitting covariance matrix. The data 
are not deep enough to provide good constraints on a, 
(j)* and M* simultaneously, and therefore the faint-end 
slopes of the LFs are assumed to be fixed at a = -0.8. 
This value is similar to the a = -0.84 ± 0.08 measured 
in the K-band for clusters at z ^ 0.3 by Muzzin et al. 
(2007a) as well as the value measured in the K-band for 
local clusters (a = -0.84 ± 0.02) by Lin et al. (2004). Al- 
though assuming a fixed value of a precludes measuring 
any evolution of the faint-end slope of the LFs with red- 
shift, it removes the strong correlation between M* and a 
in the fitting and, provided the evolution in a is modest, 
it is the best way to measure the luminosity evolution of 
the cluster galaxies via the evolution of M* . The fitted 
values of M* and the Icr errors are listed in the upper 
left of the panels in Figure 12. 

We plot the evolution of M* at 3.6/im as a function 
of redshift in Figure 14 as filled circles. Figure 14 also 
shows the predicted evolution of M* for single-burst mod- 
els with Zf = 1.0, 1.5, 2.0, 2.8, and 5.0. These models 
are normalized to M* = -24.02 at z = in the K-band, 
the result obtained by Lin et al. (2004) for 93 local clus- 
ters. This corresponds to a normalization of M* = -24.32 
at 3.6/tm, assuming a K-3.6/tm color from the z/ = 2.8 
passive evolution model. The FLS values of M* are con- 
sistent with most of these models, except the Zf = 1.0 
model, for which they are clearly too faint. Therefore, 



13 





Fig. 5. — The 3.6/jm image of the FLS with the positions of clusters superposed. The blue, green, and red circles denote clusters with 0.1 
< z < 0.4, 0.4 < z < 0.8, and z > 0.8 respectively. Large circles represent clusters with M200 > 3 x 10^* Mq and small circles represent 
clusters with M200 > 3 x lO" Mq. The size of the circles is arbitrarily chosen for clarity and is not related to the projected size of R200 
for the clusters. 



similar to the majority of previous studies we can con- 
clude that the bulk of the stellar mass in bright clus- 
ter galaxies is consistent with having been both formed 
and assembled at 2; > 1.5 and has passively evolved since 
then. As a comparison, the values measured at 3.6/im by 
De Propris et al. (2007) and Andreon (2006) are over- 
plotted as open squares and open diamonds respectively. 
These values are from spectroscopically confirmed sam- 
ples of ~ 40 clusters (the majority of which are X-ray 
detected clusters) and both agree well with the FLS val- 
ues demonstrating that passive evolution appears to be 
the ubiquitous conclusion regardless of cluster sample. 

Similar to the 3.6/xm LFs, the 4.5/zm LFs can be fit 
using a Schechter function; however, we do not perform 
fitting of the 4.5/im LFs. Instead, as a demonstration of 
the technique presented in §5.2.3 and §5.2.4, we use the 
measured 3.6/im LFs to predict the 4.5/im LFs. Unlike 
colors from the redder IRAC channels, the 3.6/im-4.5/Ltm 
colors of galaxies are nearly identical for most spectral- 



types over the redshift range 0.1 < z < 1.5. As a conse- 
quence, the 4.5/im LFs can be predicted from the 3.6/im 
LFs using the 3.6/im-4.5/im colors from almost any stel- 
lar population model. For simplicity, we use the passive 
evolution model to predict the 4.5/im LFs. The inferred 
4.5/im LFs are overplotted as solid lines in Figure 13. 
The predicted 4.5/im LFs are consistent with the mea- 
sured ones and this demonstrates that the 3.6/im LFs 
combined with simple models for the color evolution of 
galaxies can predict the LFs in other bandpasses. Fur- 
thermore, the self-consistency between the 3.6/im and 
4.5/im LFs at z = 0.15, where the 3.3/im PAH would 
contaminate the 3.6/im band, and ai z — 0.33 where it 
would contaminate the 4.5/im band suggests that the pri- 
mary source of the emission in these bandpasses at 2: < 
0.5 is stellar. 



5.2. The 5.8fim and S.Ofim Luminosity Functions 
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Fig. 6.— Multi-wavelength images of FLS J172449+5921.3 at 
Zspec = 0.252 (Cluster #10 from Table 1). The top left, top right, 
and bottom left panels arc the R-band, IRAC 3.6/im, and IRAC 
8.0/^m respectively. In each image the ficld-of-view is 1 Mpc across 
at the redshift of the cluster. The solid histogram in the bottom 
right panel shows the color distribution of galaxies with M < M* 
in the same field. The dashed histogram is the background distri- 
bution in the same aperture and the error bars show the average 
variance in the background. The arrow marks the color of the red 
sequence from the color-redshift models. The cluster red sequence 
is clearly detected at many sigma above the background. 




7.— Same as for Figure 6, but for FLS J172122-f-5922.7 at 
= 0.53 (Cluster #38 from Table 1). 



Unlike the 3.6/xm and 4.5/zm bandpasses where the 
luminosity of galaxies is dominated by emission from 
low mass stars, the luminosity of galaxies at 5.8/Lim and 
8.0/im comes from several sources. It can have contri- 
butions from warm dust continuum, PAH emission, and 
low mass stars. In particular, if warm dust (heated by 
intense star formation or an AGN) or PAH emission is 
present, it typically dominates the luminosity at these 
wavelengths. Therefore, the 5.8/im and 8.0/xm LPs can 
be useful probes of the amount of dusty star formation 
and AGN activity in clusters if the contribution from 




Fig. 8.— Same as for Figure 6, but for FLS J171648-I-5838.6 at 
Zphot = 0.56 (Cluster #44 from Table 1). 




Fig. 9.— Same as for Figmc 6, but for FLS J171420-t-6005.5 at 
Zphot = 0.63 (Cluster #50 from Table 1). 

stellar emission is properly accounted for. 

The main challenge in modeling the LPs at these wave- 
lengths is that a massive, dust-free early type galaxy 
produces relatively the same flux at 5.8/xm and 8.0/xm 
from pure stellar emission as a much lower mass star- 
burst galaxy or AGN produces from PAH emission or 
warm dust continuum. Determining the relative abun- 
dance of each of these populations in a LP is more chal- 
lenging for a statistically defined sample such as this clus- 
ter sample where individual galaxies are not identified 
as field/cluster or star forming/non-star forming. De- 
spite this challenge, we showed in §5.1 that the 3.6/im 
LPs can be used as a diagnostic of the average stellar 
emission from the cluster galaxies and that with a model 
for galaxy colors they can predict the 4.5/um LPs ex- 
tremely well. The 3.6/im-5.8/im and 3.6/im-8.0/im colors 
of different spectral-types vary significantly more than 
the 3.6/xm-4.5/im colors; however, if these colors are mod- 
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Fig. 10.— Same as for Figure 6, but for FLS J172013+5845.4 at 
Zphot = 0.69 (Cluster #56 from Table 1). 




Fig 11 — Same ds t(.i ligmc () but for FLS J171508+5845.4 at 
^phot = 0-75 (Cluster #68 from Table 1). 

clcd correctly the same technique can be used to model 
the LFs in the 5.8/xm and 8.0/xm bandpasses and pro- 
vide constraints on the number and type of star forming 
galaxies in clusters. 

Put another way, the 3.6/um LF provides effectively 
a "stellar mass budget" for predicting the 5.8/im and 
8.0/Ltm LFs. Subtracting this stellar mass budget at 
5.8/Ltm and 8.0/xm leaves an excess which can be mod- 
eled with different populations of star forming galaxies 
or AGN. Unfortunately, such models are unlikely to be 
completely unique in the sense that there will be a de- 
generacy between the fraction of star forming galaxies or 
AGN, and the intensity of the star formation or AGN ac- 
tivity within those galaxies; however, we will show that 
using only rough empircal constraints on the fraction of 
star forming/non-star forming galaxies in clusters places 
interesting constraints on the intensity of star formation 
in cluster galaxies, and the relative percentages of "reg- 



ular" star forming galaxies and dusty starbursts. 

5.2.1. Measuring the 5.8^m and S-Ojim LFs 

Before models of the cluster population are made 
we measure the 5.8/um and 8.0^m LFs using the same 
stacking and background subtraction methods as for the 
3.6/xm and 4.5/xm LFs. The LFs are plotted in Fig- 
ures 15 and 16 in the same redshift bins as the 3.6/im 
and 4.5^m LFs. IRAC is significantly less sensitive at 
5.8/im and 8.0/zm than at 3.6/im and 4.5/im and there- 
fore these LFs are much shallower. Only the bright end 
of the LF (roughly M < M* , assuming a dust-free, pure 
stellar emission early type model) can be measured with 
these data; however, this shallow depth is still sufficient 
to be a good diagnostic of the presence of luminous dusty 
starbursts. For example, at 0.1 < z < 0.4, an M82-type 
starburst is roughly 3 magnitudes brighter at 8.0/xm than 
an early type model (e.g., Huang et al. 2007a; Wilson et 
al. 2007, see also §6.3) and therefore, even a galaxy with 
M ~ M* -h 3 from the 3.6/im LF would be detected at 
8.0/xm if undergoing an M82-like dusty starburst. 

5.2.2. Contamination from AGN 

In order to draw conclusions from models of the MIR 
cluster LFs, it is important to have some constraints on 

the fraction of cluster MIR sources that are AGN and 
the fraction that are star forming galaxies. The fraction 
of galaxies in clusters at z < 0.6 identified as AGN based 
on their optical spectra in clusters is low (< 2%, e.g., 
Dressier et al. 1985, Dressier et al. 1999); whereas the 
fraction of star forming galaxies can be quite large (5 - 
80%, e.g.. Butcher & Oemlcr 1984; Dressier et al. 1999; 
Elhngson et al. 2001; Poggianti et al. 2006). There- 
fore, it might be expected that star forming galaxies will 
dominate the overall number of cluster MIR sources. It is 
possible that the AGN fraction in clusters may have been 
underestimated because some cluster AGN are missed by 
optical selection. X-ray observations of moderate redshift 
clusters have found an additional population of cluster X- 
ray AGN that do not have AGN-like optical spectra (e.g.. 
Martini et al. 2006; 2007, Eastman et al. 2007). Mar- 
tini et al. (2007) showed that this population is roughly 
as large the optical AGN population, making the overall 
AGN fraction 5% for cluster galaxies at z ~ 0.2 with 
moderate luminosity AGNs (broad-band X-ray luminosi- 
ties Lx > 10'*"'^ erg s~^), but only ^ 1% for those with 
bright AGN [hx > 10''^ erg s^^). If the analysis is re- 
stricted to galaxies with hard X-ray luminosities > 10^^ 
erg s~^, then the fraction is about an order of magnitude 
lower (0.1%, Eastman et al. 2007). 

Although these studies suggest the AGN fraction in 
clusters is low, particularly for bright AGN, it is unclear 
how many of the optical and X-ray selected cluster AGN 
will have detectable MIR emission, and what fraction 
of the cluster MIR population they comprise. Previous 
MIR studies of clusters have detected only a few AGN 
in spectroscopic samples of ^ 30-80 cluster MIR sources 
(e.g., Due ct al. 2002; Coia et al. 2005; Marcillac et al. 
2007; Bai et al. 2007) suggesting that >90% of cluster 
galaxies detected in the MIR are star forming galaxies. 
One way to estimate the fraction of cluster MIR-bright 
AGN is to use the IRAC and MIPS color-color diagrams 
suggested by Lacy et al. (2004) and Stern et al. (2005). 
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Fig. 12. — The 3.6/Lun LFs of clusters in the FLS. The solid line shows the best-fit Schechter function assuming a = -0.8. The redshift, 
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Fig. 14. — Evolution in M* from the 3.6/ini LFs as a function 
of redshift. The long dashed, dash-dotted, dash-dot, solid, and 
dashed lines show models where the stars form in a single burst at 
z = 1.0, 1.5, 2.0, 2.8, and 5.0 respeetively. The filled circles are 
the FLS clusters and the open diamonds and open squares are the 
M* values from the Andreon (2006) and De Propris et al. (2007) 
cluster samples. 

Although these simple color cuts fail to identify com- 
plete samples of AGN because they only identify those 
that have red power-law slopes in the MIR (e.g., Carda- 
mone et al. 2008); these are precisely the type of AGN 
that will be included in the 5.8/im and 8.0/Ltm LFs and 
therefore the color cuts should provide a reasonable esti- 
mate of the contamination of those LFs from AGN. 

In the left panels of Figure 17 we plot the IRAC colors 
of all galaxies brighter than the 50% completeness limits 
using the color spaces suggested by Stern et al. (2005) 
(top) and Lacy et al. (2004) (bottom). The dashed lines 
in each panel represent the portion of color space used to 
select AGN in the MIR by these authors. FLS galaxies 
which satisfy the color criteria are plotted as grey circles. 
The right panels of Figure 17 shows the same plot for all 
galaxies with R < R200 for clusters at < 0.7 in the FLS 
(59 clusters). 

The entire FLS (left panels) can be used to estimate 
the surface density of MIR selected AGN in these color 
spaces. Subtracting this background from the cluster 
fields w(^ find and excess of 26 ± 22 galaxies using the 
Stern ct al. (2005) color cut and and excess of 30 ± 30 
galaxies using the Lacy ct al. (2004) color cut. Summing 
the background subtracted 3.6/xm LFs to the same limit 
implies that there are 2466 total cluster galaxies in these 
59 clusters, and that the overall fraction of cluster galax- 
ies that candidate MIR-bright AGN (to our 3.6^m detec- 
tion limit) is lli[%, where all error bars have been calcu- 
lated using Poisson statistics. Integrating the 5.8/Ltm and 
8.0/xm LFs shows there are 869 and 959 cluster galaxies 
detected in these bands and that the fraction of cluster 
sources detected in the MIR that are candidate AGN is 
~ 31^%. 

Although this crude estimate is almost certainly an 
incomplete census of the total fraction of AGN in clus- 
ters, it is remarkably similar to the AGN fractions mea- 
sured with optical spectroscopy or by X-ray selection and 
is consistent with the fraction of spectroscopically con- 
firmed MIR-bright AGN seen in previous cluster MIR 
studies. Based on the low estimated AGN fraction, and 
for the sake of simplicity in interpretation, we do not 
model an AGN component in the 5.8/im and 8.0/xm LFs 
in this analysis. We do note that the X-ray, spectroscopic 



and MIR selection do show clearly that the fraction of 
MIR cluster sources that are AGN is not zero, and there- 
fore some of the sources in the 5.8/um and 8.0/xm LFs will 
certainly be AGN. 

5.2.3. Modeling the 5. Sum Luminosity Function 

The simplest fiducial model that can be made for the 
MIR cluster galaxy population is to assume that the 
bright end of the LF is dominated by passive, dust-free, 
early type galaxies (i.e., the emission at 5.8/Ltm and 8.0/im 
is completely stellar). Although such a model is unreal- 
istic, it provides a baseline for predicting the amount of 
emission in the MIR from stellar emission, and any ex- 
cess beyond this model is likely to be from dusty star 
formation in the cluster population. Assuming such a 
model, the 5.8/um LFs can be inferred from the 3.6^m 
LFs using the 3.6/im-5.8/im colors from the Bruzual & 
Chariot passive evolution model. These predicted 5.8/im 
LFs arc ovcrplottcd on the LFs in Figure 15 as the solid 
red lines (Figure 15 also has additional models overplot- 
ted which are introduced in §5.2.4). 

Qualitatively, the 3.6/im LFs and the passive evolu- 
tion model predict the 5.8/zm LFs reasonably well at all 
redshifts. This is perhaps not surprising because due to 
k-corrections, 5.8/tim is only sensitive to emission from 
warm dust or PAHs in star forming galaxies at z < 0.3 
(see §6.3). For galaxies at higher redshift. 5.8/im probes 
rest-frame wavelengths which, similar to the 3.6/im LFs, 
are dominated by stellar emission. As a result, any dusty 
star forming cluster galaxies would only be visible as a 
notable excess in the predicted 5.8/im LFs at 2; < 0.3. 
No such excess is seen: however, the fraction of blue 
star forming galaxies in clusters evolves rapidly (i.e., the 
Butcher-Oemler Effect) and clusters at z < 0.3 typically 
have low blue fractions and relatively few star forming 
galaxies (e.g., EUingson et al. 2001; Balogh et al. 1999; 
Margonincr ct al. 2001). This result confirms that the 
fraction of star forming galaxies in clusters at z < 0.3 low, 
and, that furthermore there is no significant additional 
population of MIR luminous dusty star forming galax- 
ies in clusters at these redshifts that are missing from 
optically-selected spectroscopic or photometric studies. 

5.2.4. Modeling the S.Ofim Luminosity Function 

Unlike the 5.8/um LFs, the cluster 8.0/im LFs are not 
consistent with the passive evolution model predictions 
from the 3.6/im LFs illustrated by the solid red lines plot- 
ted in Figure 16. This model clearly underprcdicts the 
number of galaxies in the 8.0/im LFs at all redshifts. In 
order to construct a more useful model for the 8.0/im 
LF that includes the cluster star forming population, we 
use the 3.6/im-8.0/im colors for different types of star 
forming galaxies from J. Huang et al. (2008, in prepa- 
ration). These authors have empirically extended the 
color /redshift models of Goleman et al. (1980) to 10/im 
using local galaxies with ISO spectroscopy. Some ex- 
amples of the colors from these models are presented in 
Wilson et al. (2007). 

Given the large nunibca' of permutations possible in 
the types of star forming galaxies, we are interested in as 
simple a model as possible which will allow for a straight- 
forward interpretation of the data. For this analysis 
we divide the cluster star forming population into two 
populations: "regular" star forming cluster spirals, and 
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Fig. 15. — The 5.8/im LFs of clusters in the FLS. The solid red line shows the 5.8/im predicted from the 3.6/im LF assuming all 
galaxies have the colors of the passive evolution model. The dashed green lines and dotted blue hnes are the regular+quiescent and 
starburst+regular+quiescent models described in §5.2.4 respectively; however, 5.8/im is not sensitive to PAH emission or warm dust at z > 
0.3 and therefore these models are not notably different from the passive evolution model. 



dusty starburst galaxies. Huang et al. (2008) have mod- 
els for both Sbc and Sod galaxies; however, the 3.6/im- 
8.0/ini colors of these models are indistinguishable, and 
therefore we adopt their Sbc galaxy as the model for 
a "regular" star forming cluster spiral. Huang et al. 
also have colors for several "canonical" dusty starburst 
galaxies such as M82, Arp220, and NGC 1068. M82 is 
a moderate-strength dusty starburst, has no AGN com- 
ponent, and is classified as a luminous infrared galaxy 
(LIRG). By contrast, Arp220 and NGC 1068 are power- 
ful dusty starbursts with AGN components. The IR lu- 
minosity of Arp220 is dominated by star formation from 
a major merger, while the IR luminosity of NGC 1068 is 
dominated by a powerful AGN (although both galaxies 
have AGN and starburst components). Both are classi- 
fied as ultra- luminous infrared galaxies (ULIRGs). Given 
that the majority of distant clusters studied thus far in 
the MIR have shown a significant population of LIRGs 
but no population of ULIRGs (e.g., Coia et al. 2005; 
Geach et al. 2006; Marcillac et al. 2007), we assume 
that any cluster dusty starbursts will have colors similar 
to M82, rather than Arp220 or NGC 1068. In general, 
replacing M82 as the model for cluster dusty starbursts 
with cither Arp220 or NGC 1068 requires a smaller frac- 
tion of dusty starbursts since they are more luminous. 

In order to ascertain the dominant mode of star forma- 
tion present in the cluster population we can construct 
simple models for the 8.0/im LFs from the 3.6/im LFs 
using various combinations of these populations. The 
purpose of the models is not to perfectly reproduce the 
cluster 8.0/im LFs (this requires a much more detailed 



knowledge of the populations in each cluster than can be 
obtained by statistical background subtraction), but to 
demonstrate how the 8.0/im LFs should appear given dif- 
ferent proportions of these populations and thereby esti- 
mate the importance of each's contribution to the 8.0/im 
LFs. Hereafter we refer to the Sbc model as "regular", 
the M82 model as "dusty starburst" , and the Bruzual & 
Chariot passive evolution model as "quiescent" . 

Beyond assuming that all cluster galaxies are quies- 
cent, which clearly underpredicts the 8.0/im LFs, the 
next most simple model that can be made is to as- 
sume some fraction of the cluster galaxies are "regu- 
lar" star forming galaxies (hereafter we refer to this 
model as regular-|-quiescent). In order to make such a 
model we require an approximation of the relative pro- 
portions of star forming and quiescent galaxies in clus- 
ters as a function of redshift and luminosity. The best 
spectroscopically-classified data at these redshifts comes 
from the MORPHS (Dressier et al. 1999; Poggianti et 
al. 1999) and CNOCl (Balogh et al. 1999; EUingson et 
al. 2001) projects. Unfortunately, the number of cluster 
spectra per dz is relatively small in these samples and 
they cover only a modest range in redshift (0.2 < z < 
0.5) and depth in terms of the cluster M*. 

Although spectroscopic classification would be the 
most reliable, the lack of data motivates the use of clus- 
ter blue fractions (ff,) as a function of redshift as a model 
for the relative fractions of star forming/non-star form- 
ing galaxies. Blue fractions for reasonably large samples 
of clusters at different redshifts have been calculated and 
it is fairly straightforward to measure them as a function 
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Fig. 16. — The 8.0/jm LFs of clusters in the FLS. The soUd red line, dashed green line, and dotted blue line are the S.O^m LFs predicted 
using the 3.6/^m LF and the quiescent, regular+quiescent, and starburst+regular+quiescent models described in §5.2.4. At lower redshift 
(z < 0.4) the LFs are most similar to the predictions from the rcgular+quiesccnt model whereas at higher redshift (z > 0.4) the LFs are 
better described by the starburst+regular+quiescent model. 

some variation in as a function of magnitude, all of the 
model LFs consistently overprcdict the number of bright 
galaxies in the 8.0/im LFs, and underpredict the num- 
ber of faint ones. In order to estimate the variation of f^ 
as a function of magnitude we use the spectrally-typed 
LFs of Muzzin et al. (2007a). They measured the K- 
band LF for cluster galaxies defined spectroscopically as 
either star forming or quiescent. Comparing those LFs 
(their Figure 13) and assuming all star forming galaxies 
are blue, and all quiescent galaxies are red, results in f^ 
values of 0.19, 0.35, and 0.52 for galaxies with M < M*, 
M* < M < M* 1, and M* +1 < M < M* + 2 re- 
spectively in clusters at z ~ 0.3. Comparing these values 
shows that % is 1.8 times larger at M* < M < M* + 1 
than at M < M*, and is 2.7 times larger at M* +1 < 
M < M* + 2 than at M < M*. We therefore adopt an 
fb that varies with magnitude with the following condi- 
tions: For galaxies with M < M* in the 3.6/xm LF we use 
the ffc values from Table 2. For galaxies with M* < M < 
M* + 1, we assume that fb is twice as large as the values 
in Table 2, and for galaxies with M* + 1 < M < M* + 
2 we assume that fh is three times as large as the values 
in Table 2. In cases where this causes ffc > 1.0, it is set 
equal to 1.0. 

Combining the f{, as a function of redshift and mag- 
nitude with the 3.6/L(m LFs assuming all "blue" galaxies 
have the color of the Huang et al. Sbc galaxies and all 
"red" galaxies have the color of the passive evolution 
model results in the models that are plotted as green 
dashed lines in Figures 15 and 16. Comparing the data to 
these models shows that this simple model using only reg- 



of magnitude within these clusters. In particular, using 
ffc as an estimate of the star forming fraction should pre- 
dict the number of blue star forming galaxies (i.e., those 
with colors similar to the Sbc model). If a population 
of red, dust-obscured starburst galaxies exists in clusters 
they should be evident in the 8.0/zm LFs as an excess of 
galaxies beyond the regular+quiescent model. 

For ff, as a function of redshift we use the data of EUing- 
son et al. (2001) from the CNOCl clusters which span 
the redshift range z — 0.2 to z = 0.4, and for clusters 
at z > 0.4 we use the data on RCS-1 clusters from Loh 
et al. (2008). Rough ff, values for both these samples 
were recomputed using only galaxies with M < M* (D. 
Gilbank private communication), because this matches 
the depth of the 5.8/im and 8.0/.tm LFs. These f;, values 
as a function of redshift are listed in Table 2. 

The scatter in cluster f;, values at a given redshift is 
large, and therefore different studies find different mean 
values depending on sample. The values we have adopted 
are consistent with the majority of work in the field (e.g.. 
Butcher & Oemler 1984; Small et al. 1998; Margoniner 
et al. 2001; Andreon et al. 2004), although we have mea- 
sured them using a brighter luminosity cut. Of course, 
the best way to infer the f^ of the FLS clusters would be 
to measure it from the clusters themselves; however, we 
do not have the proper filter coverage at 2; < 0.5 to make 
this measurement properly nor a large enough sample to 
make a measurement that would be statistically different 
from the adopted values. 

The cluster f^ is also a function of limiting magnitude 
(e.g., EUingson et al. 2001), and without incorporating 
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Fig. 17. — Top Left Panel: Color-color plot of all galaxies in the FLS (small dots). The dashed lines denote the region used to select 
AGN by Stern et al. (2005). Bottom Left Panel: Same as top left, but for the Lacy ct al. (2004) color space. Right Panels: Color-color 
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at 3.6/im are AGN and that ~ 3^3% of cluster galaxies detected at 5.8/im and 8.0/im are AGN. 



ular-|-quicsccnt galaxies predicts the cluster 8.0/ini LFs 
fairly well. In particular, the z = 0.15, 0.25 and 0.33 LFs 
are well described by this model. For the higher redshift 
LFs this model is clearly better than the purely quiescent 
model; however, it still does not account for the entire 
8.0/um population. 

Most importantly, the regular-|-quiescent model shows 
that out to 2; ~ 0.65, where S.O^m still probes rest-frame 
dust emission, there is no significant population of bright 
(M < M*) galaxies in clusters that cannot reasonably be 
accounted for by "regular" star forming cluster spirals. 
This is significant because it suggests that whatever pro- 
cesses responsible for transforming the morphology and 
spectral-type of bright cluster galaxies over the same red- 
shift range do not involve an ultra-luminous dusty star- 
burst phase such as those caused by major mergers of 
gas-rich galaxies (i.e., "wet" mergers). We note that 
there appears to be an overdensity of very bright galax- 
ies in the z = 0.82 LF that cannot be accounted for by 
the regular+quiescent model and this suggests the pos- 
sibility of an onset of luminous starbursts (possibly from 
mergers) or AGN activity in bright galaxies at higher 
redshift. 

Although the regular -l-quiescent model predicts the 

bright end of the 8.0/im LFs well at all rcdshifts, and 
the entire 8/um LF at lower redshift, it fails to account 



for all of the LFs. In particular, this model seems to 
underpredict the number of fainter galaxies in the 8.0/im 
LFs for clusters at z > 0.4. This suggests a third com- 
ponent to the cluster 8.0/im population, possibly a red, 
dusty starburst population which is not accounted for 
by the cluster ff,. Such a population was suggested by 
Wolf et al. (2005) who found that the SEDs of roughly 
30% of the red sequence galaxies in the Abell 901/902 
supercluster {z = 0.17) were better described by dusty 
templates rather than a dust-free, old stellar population. 
In order to explore this possibility, we construct a new 
model with the same values of ff, as a function of mag- 
nitude and redshift as for the regular -l-quiescent model, 
but this time we assume that some of the red quiescent 
galaxies arc instead M82-like dusty starbursts. M82 has 
optical-IR colors that are similar to quiescent galaxies 
(see Huang ct al. 2008 and §6.3) so it is reasonable to 
assume that any M82-like dusty starbursts would be part 
of the population of red cluster galaxies rather than the 
blue cluster galaxies. 

If we assume that the dusty starburst population is a 
constant fraction of the red cluster galaxies, this would 
result in a varying ratio of dusty starburst to regular star 
forming galaxies in clusters as a function of redshift. In 
particular, clusters at low redshift will have the highest 
fraction of dusty starburst galaxies (because the fb is low 
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and the red fraction is high) . The LFs above have already 
suggested that there is no need for a dusty starburst pop- 
ulation at low redshift, so modeling the dusty starbursts 
as a fixed fraction of the red galaxies seems inappropri- 
ate. Instead, a better way to model the population is to 
assume that the cluster ff, is a tracer of the total star for- 
mation in the cluster and that ratio of dusty starburst to 
regular star forming galaxies is a constant. Given this as- 
sumption we can predict the fraction of dusty starbursts 
directly from the cluster 4. This fraction of dusty star- 
bursts is then removed from the fraction of red quiescent 
galaxies and a model for the LFs can be made. Hereafter 
we refer to this model as starburst -|-regular-|-quiescent. 
The fractions of the cluster galaxy populations in terms 
of ib are defined using the equations, 

fdsb = Ib X fdsb/reg, (4) 
fq = ^- fb- fdsb, (5) 

where idsb is the fraction of dusty starburst galaxies, 
^dsb/reg is the assumed ratio of dusty starburst to regular 
star forming galaxies, and fg is the fraction of quiescent 
galaxies. In cases where fdsb + f& > 1 we set f^sh = 1 - fb 
and iq = 0. 

As of yet there are no good observational constraints 
on the parameter fdsb/reg- Therefore, as a first-order fidu- 
cial value we assume that fdsb/reg = 0.5. In general, we 
find that allowing a range of values between 0.3 - 1.0 pro- 
vides models that are fairly similar. More importantly, 
the differences in models that use igb/reg between 0.3 - 
1.0 are much smaller than the difference between any of 
those models and the regular -(-quiescent model. There- 
fore, the interpretation of the data using these models 
will not depend strongly on the assumed value of fst/reg- 
The star burst -|-regular-|-quiescent model with isb/reg = 
0.5 is overplotted on Figures 15 and 16 as the dotted 
blue line. 

This starburst -l-regular-l-quiescent model over-predicts 

the number of bright galaxies in the z < 0.4 8.0/im LFs 
but it is better at describing the LFs a,t z > 0.4 than 
the regular -l-quiescent or purely quiescent models. This 
suggests that there is a population of dusty starbursts 
in clusters at 2; > 0.4 that does not exist at 2; < 0.4, 
and that these starbursts arc consistent with being of an 
M82-type. We discuss this in more detail in §6.1. 

6. DISCUSSION 

6.1. Evidence for a Change in Star Formation 
Properties of Cluster Galaxies? 

In order to better illustrate the differences in the model 
populations described above, we subtract the quiescent 
model from the 5.8/im and 8.0/im LFs between 0.15 
< z < 0.65 and plot the residuals in Figures 18 and 
19. The residuals from the quiescent H-regular model and 
starburst -(-regular -hquiescent models from §5.2.4 are also 
plotted in Figures 18 and 19. The solid vertical lines in 
the plots represent the magnitude of M* inferred from the 
3.6/im LF assuming the passive evolution model, and give 
some indication of the depth of the LFs. If we compare 
the data to the models and take the results at face value, 
it suggests that the intensity of star formation in clusters 
is evolving with redshift and that it can be classified into 
three types. The first type of star formation is "weak" 
and best describes the lowest redshift clusters {z < 0.15) 



which are consistent with the colors of an almost ex- 
clusively quiescent population in all IRAC bandpasscs. 
This result is consistent with numerous studies of nearby 
clusters using spectroscopy which show few star forming 
galaxies (e.g.. Dressier et al. 1985, Popesso et al. 2007; 
Christlein & Zabludoff 2005; Rines et al. 2005). 

Between 0.2 < z < 0.5 the 8.0/im LFs arc no longer 
well-described by the purely quiescent model and the reg- 
ular -(-quiescent model is the best model. This shows that 
the majority of star formation in clusters at this epoch 
is primarily relegated to galaxies that have MIR colors 
similar to local late- type star forming galaxies (i.e., the 
Sbc model). This has direct implications for the SFRs of 
these galaxies because Wu et al. (2005) showed that the 
dust-obscured SFR of galaxies is proportional to their 
8.0/im flux. Although other authors have demonstrated 
that there arc caveats when rising the 8.0/im flux to infer 
SFRs (i.e., the scatter can be as high as a factor of 20- 
30, Dale et al. 2005), this still implies that the average 
SFR or the average SFR per unit stellar mass (the av- 
erage specific star formation rate, SSFR) of star forming 
cluster galaxies at 0.2 < z < 0.5 is similar to those in 
the local universe (because they have 3.6/im-8.0/im col- 
ors similar to local Sbc galaxies). This second mode of 
star formation in clusters is roughly what would be con- 
sidered "regular" star formation for galaxies in the local 
universe. 

At 2: > 0.5 the starburst+rcgular+quiesccnt model be- 
comes the best description of the LFs. Again, assuming 
that 8.0/xm flux is an indicator of SFR, the M82 starburst 
model is approximately a factor 2.5 brighter at S.O/zm 
than the regular Sbc model for the same 3.6/im flux. 
Given that our model suggests that regular star form- 
ing galaxies make up ^ 30-40% of the cluster population 
at this redshift and M82 galaxies make up ~ 15-20%, 
this implies that not only is the abundance of star form- 
ing galaxies in clusters higher at higher redshift (i.e., the 
Butcher-Oemler Effect), but also the average SSFR of 
cluster galaxies is approximately a factor of 1.5 higher at 
z > 0.5 than it is at 2 < 0.5. This increase in SSFR sug- 
gests a third mode of star formation in cluster galaxies 
that could be considered a "burst" mode, at least relative 
to local star formation rates. Interestingly, this increase 
in the SSFR of cluster galaxies at higher redshift is con- 
sistent with field studies of the universal star formation 
density (p*) which show an increase of roughly a factor of 
2-5 between z = 0.2 and z = 0.5 (e.g., Lilly et al. 1996; 
Madau et al. 1996; Wilson et al. 2002; Scliiminovich 
et al. 2005; Lc Floc'h et al. 2005). It suggests that the 
increasing fraction of dusty starbursts in the cluster pop- 
ulation could be interpreted as the result of an increase in 
the universal SSFR of galaxies with redshift and the con- 
stant accretion of these galaxies into clusters and is not 
necessarily because starbursts arc triggered by the clus- 
ter environment. Furthermore, these galaxies might only 
be considered "starbursts" relative to the mean SSFR 
locally, whereas at higher redshift their higher SSFR is 
simply typical of galaxies at that redshift. We compare 
the cluster 5.8/im and 8.0/im LFs to the field LFs in §6.2 
and discuss this further in that section. 

It is interesting that the cluster star forming popula- 
tion transitions from being best described by regular star 
forming galaxies to regular and dusty starburst galax- 
ies around a redshift of 2; 0.4. This is notable be- 
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cause of the discrepant abundances of k+a and a+k post- 
starburst galaxies found in clusters by the MORPHS 
(Dressier at al. 1999) and CNOCl (Balogh et al. 1999) 
projects. Dressier et al. (1999) found that approximately 
18% of cluster galaxy spectra could be classified as k+a 
galaxies based on the equivalent width of the H(5 line, 
whereas Balogh et al. (1999) found that only 2% of the 
cluster population could be classified this way. These re- 
sults obviously lead to very different interpretations of 
the role of starbursts in the evolution of cluster galaxies. 
In particular, Dressier et al. found that the number of 
k+a galaxies was an order of magnitude higher in clus- 
ters than the coeval field, suggesting a cluster-related 
process to the creation of these galaxies, while Balogh et 
al. found roughly equal numbers, suggesting no environ- 
mental role. 

Although both Dressier et al. (2004) and Balogh et 
al. (1999) have pointed out that the different methods 
of data analysis may be partly responsible for such dis- 
crepant results, this study suggests that the slightly dif- 
ferent redshift range of the MORPHS and CNOCl sam- 
ple may also play some role. Excluding the two highest 
redshift clusters in the CNOCl sample (MS 0451-03 and 
MS 0016+16, both at z ~ 0.55) the mean redshift of the 
other 14/16 (88%) clusters in the sample is z = 0.28. By 
contrast, the mean redshift of the MORPHS sample is 
z = 0.46. Our 8.0/im cluster LFs seem to indicate that 
z ^ 0.4 represents a transition redshift above which the 
dominant mode of star formation in clusters is better de- 
scribed as starburst, as opposed to regular. Given that 
once star formation ceases, the typical lifetime of the A 
star component of a starburst galaxy's spectrum is ^ 1.5 
Gyr, and that the lookback time between z — 0.46 and 
z = 0.28 is also 1.5 Gyr, it is possible that both dusty 
starbursts, and k+a galaxies that are in clusters at ^; = 
0.46 may have evolved to quiescent "k"-type galaxies by 
z ^ 0.28, provided that the dusty star formation is im- 
mediately truncated. This would be consistent with the 
change in the 8.0/im LFs around this redshift and may 
explain why the MORPHS and CNOCl samples show 
different abundances of post-star burst galaxies. Further- 
more, 1.5 Gyr prior to z = 0.46 is z ^ 0.65. Our z = 0.65 
cluster LF has the largest abundance of dusty starburst 
galaxies, and if a significant fraction of these had their 
star formation truncated, these would be logical progen- 
itors to the large population of k+a galaxies seen at z = 
0.46 by Dressier et al. (1999). 

Our results, which show an increase in the strength of 
the dominant mode of star formation in cluster galax- 
ies (from weak to normal to starburst), as well as an 
overall increase in the abundance of dusty star form- 
ing galaxies are also consistent with MIR observations 
of other clusters at these redshift ranges. In particular, 
Coia et al. (2005), Geach et al. (2006), Marcillac et al. 
(2007), and Bai et al. (2007) have all shown that clusters 
at higher redshifts have significantly more MIR sources 
than clusters at lower-redshift, and that these sources are 
typically brighter than the sources in lower-redshift clus- 
ters. Taken at face value, our results and their results 
show the equivalent of a Butcher- Oemler Effect in the 
MIR where both the fraction, and SSFR of star forming 
galax;ies is increasing with increasing redshift. Whether 
this increase is caused by the increase in the universal 
SFR with redshift, and the constant infall of such galax- 



ies into the cluster environment, or by the triggering of 
starbursts by the high-redshift cluster environment is still 
uncertain. We investigate this point further in §6.2 by 
comparing the cluster and field IRAC LFs. 

6.2. 7s the Cluster Population Different From the Field 
Population? 

The most obvious way to understand if the cluster en- 
vironment is responsible for triggering starburst events is 
to directly compare the field and cluster 5.8/im or 8.0/im 
LFs and look for an excess of galaxies in the cluster LFs. 
For this comparison we use the field LFs measured by 
Babbedge et al. (2006, hereafter B06). Their LFs are de- 
termined using photometric redshifts of ~ 100 000 galax- 
ies from a 6.5 deg^ patch of the SWIRE survey. The field 
LFs are measured in 5 redshift bins, and we compare the 
cluster LFs to the three bins which overlap the redshift 
range of the clusters (0.0 < z < 0.25, 0.25 < z < 0.50, 
and 0.5 < z < 1.0). The corresponding cluster LFs used 
for comparison are the z = 0.15, z = 0.33, and z = 0.65 
LFs respectively. 

The B06 field LFs are determined using total luminosi- 
ties, not apparent magnitudes like for the cluster LFs. 
Converting the units of the cluster LFs to total lumi- 
nosities requires distance moduli and full k-corrections. 
In §5.2.4 we showed that the cluster LFs can be well- 
described using three basic populations of galaxies: qui- 
escent, regular star forming, and dusty starburst. We 
use the models of these three spectral types for the k- 
corrections. The k-corrections for the quiescent galax- 
ies are taken from the single-burst model and the k- 
corrections for the regular and dusty starburst galax- 
ies arc taken from the Huang et al. (2008) Sbc and 
M82 models respectively. Each LF is statistically k- 
corrected using the relative proportions of the galaxies 
which best described the LFs in §5.2.4. The apparent LF 
for each redshift is divided into the three components by 
the fraction of galaxies of that type and are individually 
k-corrected and shifted by the distance modulus. These 
LFs are then summed to provide the total cluster LF in 
terms of absolute luminosities in units of vLi,/Lq. 

The cluster LFs are normalized by the number of galax- 
ies per virial volume, whereas the field LFs are normal- 
ized by their actual number density per Mpc^. The clus- 
ter normalization can be put in the same units as the 
field LFs by dividing by the virial volume; however, this 
does not provide a fair comparison because clusters have 
much higher volume densities of galaxies than the field. 
The most useful way to compare the cluster and field 
LFs is on a per unit stellar mass basis. We do not have 
stellar mass functions for either the field or cluster; how- 
ever, we can again assume that the 3.6/im himinosity is 
roughly a proxy for stellar mass and renormalize the LFs 
to a common normalization so that they reproduce the 
same </>* in the Schechter function fits. The renormalized 
3.6/im, 4.5/im, 5.8/im, and 8.0/im cluster LFs are plotted 
in Figures 20, 21, 22, and 23, respectively as the filled 
red circles. The field LFs are ovcrplotted as blue squares. 

Figures 20 and 21 show that the overall shape of the 
cluster and field 3.6/im and 4.5/im LFs are similar at all 
redshifts. There is a slight, though not statistically sig- 
nificant, excess in the number of the brightest galaxies in 
the cluster LFs; however, these are likely to be giant ellip- 
tical galaxies which are common in clusters and typically 
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do not follow the distribution of the Schechter function. 
Other than the giant ellipticals, the shape of the 3.6/im 
and 4.5/ini cluster and field LFs are similar which shows 
that the distribution of galaxies as a function of stellar 
mass is nearly identical in these environments. This re- 
sult is consistent with K-band studies which have shown 
only small differences in M* (< 0.2 mag) between these 
environments (e.g., Balogh et al. 2001; Lin et al. 2004; 
Rines et al. 2004; Muzzin et al. 2007a). 

Conversely, there are significant differences in the 
5.8/im and 8.0/im LFs of the cluster and field. Both the 
5.8/im and S.O/xm LFs follow a sequence where the cluster 
LF is more abundant in MIR galaxies at z = 0.65, partic- 
ularly moderate- luminosity galaxies, and thereafter the 
abundance of MIR galaxies in clusters declines relative to 
the field with decreasing redshift. At z = 0.33, the clus- 
ter is slightly deficient in both 5.8/im and 8.0/im galaxies 
relative the field, reduced by a factor of ~ 2 for galax- 
ies with i/L^ = 5 X 10^ - 5 X 10^° Lq. At z = 0.15, 
the cluster LF is significantly depleted compared to the 
field, reduced by a factor of ~ 5 for galaxies with i^hi, 
= 5 X 10^ - 5 X 10^° Lq. This trend not only indicates 
that the environment of dusty star forming galaxies af- 
fects their evolution, but that the enironmental effects 
seem to evolve with redshift. At z = 0.15 dusty star 
forming galaxies are more frequently found in the lower 
density field environment, whereas at z = 0.65 they are 
found more frequently in the higher density cluster envi- 
ronment. 

Our results are similar to those from recent studies by 
Elbaz et al. (2007) and Cooper et al. (2008) that have 
shown that the mean star formation rate of field galax- 
ies in higher density environments increases faster than 
those in low density environments with increasing red- 
shift. This differential increase leads to a remarkable re- 
versal in the slope of the <SFR> of galaxies as a function 
of density at z ~ 1 as compared to z ~ 0. Field galax- 
ies in high density environments at z ~ 1 actually have 
higher <SFR> than those in low density environments. 
Although those studies compare <SFR> of galaxies at a 
range of densities within the field and do not use clusters 
per se, our comparison between the 8/im LFs of the clus- 
ter and field environments seem to at least qualitatively 
suggest a similar trend. 

It is not entirely obvious why starbursts should pre- 
fer the cluster environment over the field environment at 
high (z > 0.5) redshift and then reject it at lower red- 
shift (z < 0.5). We suggest that starbursts could pref- 
erentially be triggered during the initial formation and 
collapse of the cluster, and be quenched thereafter by 
the high-density environment. If this interpretation is 
correct, it is likely that the parameter most responsible 
for the change in star formation properties relative to the 
field is the degree of virialization of the clusters. 

Clusters that are unrelaxed, or in the process of col- 
lapsing, have two properties that would permit increased 
numbers of dusty starbursts. Firstly, before virializa- 
tion, the cluster gas has not yet been shock-heated to 
its maximum temperature. This hot intracluster gas has 
long been considered the primary cause for the quench- 
ing of star formation in cluster galaxies because it pre- 
vents the cooling of gas in the outer halo of a galaxy, 
thereby "strangling" star formation. Depending on the 
density/temperature threshold required for quenching. 
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Fig. 20. — Comparison between the cluster and field 3.6/im LFs 
at diff'erent redshifts. The field LFs are plotted as open blue squares 
and the cluster LFs are plotted as filled red circles. The cluster 
LFs are renormalized so that the values of ij>* from the Schechter 
function fits (§5.1) match the 0* values from the Schechter function 
fits in B06. 
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Fig. 21. — Same as Figure 19 but for the 4.5/im LFs. 




Fig. 22. — Same as Figure 19 but for the 5.8/(m LFs. 
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Fig. 23. — Same as Figure 19 but for the 8.0/im LFs. 



it is possible that starbursts that would normally be 
quenched in virialized clusters at lower redshifts may sur- 
vive longer in unvirialized clusters at high redshift. Sec- 
ondly, the velocity dispersions in unrelaxed systems are 
lower and therefore mergers and harassments should be 
more common at higher redshift (e.g., Tran et al. 2005b). 
It is plausible that this more dynamically "active" en- 
vironment preferentially triggers star formation. The 
combination of more triggered dusty starbursts through 
harassment and mergers as well as a weaker quenching 
process may be the reason for more dusty starbursts in 
clusters relative to the field at higher redshift. Once a 
cluster becomes virialized the interactions between galax- 
ies should become less frequent and the quenching of star 
formation by the hot cluster gas will be more efficient. 
In such a scenario the relative abundances of dusty star- 
bursts in clusters should decrease relative to the field. 

If our interpretation is correct we might expect differ- 
ent results from the 8.0/im LFs of X-ray selected samples 
of clusters (i.e., those which require a hot virialized clus- 
ter gas component) compared to red sequence selected 
samples, which, assuming the early type population is 
formed prior to cluster collapse, do not require that clus- 
ters are fully virialized. 

6.3. Are the Color Models Correct? 
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The main conclusions from the cluster LFs presented in 
this paper depend on interpreting color models that have 
been primarily calibrated or determined using nearby 
galaxies. If these models are not applicable at higher 
redshift then this could cause incorrect conclusions to be 
drawn from the LFs. Using the spectroscopic redshifts 
we can examine the colors of confirmed cluster galaxies 
as a function of redshift to check if the models are rea- 
sonable. 

There are 55 spectroscopic redshifts available for clus- 
ter galaxies (sec §2.4 & §2.5). Using the spectra we can 
classify these galaxies into two basic types, star form- 
ing and non-star forming. For the Hectospec, SDSS, 
and WIYN spectroscopy the best-fitting cross-correlation 
template is used for the classification. For the remaining 
galaxies the classification is made by eye-examining the 
spectra for any evidence of the [Oil] , [OIII] , or Ha emis- 
sion lines. Galaxies with any of these emission lines are 
classified as star forming, and those without are classified 
as non-star forming. Although this is a crude approach 
to classifying galaxies, we are only interested in a rough 
classification and taking a more quantitative approach, 
such as measuring EWs, is unnecessary. Furthermore, in 
all cases the cluster galaxies had spectra that were typ- 
ical of either normal star forming (several emission lines 
including [Oil] and Ha) or quiescent galaxies (strong H 
and K lines and a 4000A break), and classification was 
straightforward. There were no hybrid objects associated 
with clusters except two AGN from the Hectospec data. 

In Figure 24 we plot several of the colors of these 
galaxies as a function of redshift. Star forming galax- 
ies are plotted as purple points and non-star forming 
galaxies are plotted as red points. The Bruzual & Char- 
lot single-burst model is overplotted as the solid line, 
and the Huang et al. (2008) Sbc and M82 models are 
overplotted as the dotted and dash dotted lines, respec- 
tively. In general, the non-star forming galaxies follow 
the single-burst model well at all redshifts. There are a 
handful non-star forming galaxies which appear to have 
some excess B.O/xm emission, and this may be from either 
low-level star formation or a low-luminosity AGN. 

There are fewer star forming than non-star forming 
galaxies in the sample; however, their colors follow the 
Sbc and M82 models quite well. At 8.G/im, where the col- 
ors of the Sbc and M82 models are most different from 
the single-burst model, it is clear that galaxies with emis- 
sion lines have colors similar to those models, whereas 
those without tend to follow the single-burst model. Half 
of the star-forming galaxies in Figure 24 (8/16) come 
from our spectroscopy of FLS J172449-I-5921.3 (cluster 
#10, z = 0.252). These galaxies were selected for spec- 
troscopy because they were detected at 24/im. Interest- 
ingly, most of these galaxies (7/8) have a 3.6/im - 8.0/im 
color similar to the Sbc model, yet they show a wide 
range in R - 3.6/um color. A few have an R - 3.6/xm 
color bluer than the red sequence, typical of Sbc galaxies, 
whereas others have an R - 3.6/im color redder than the 
red sequence. This illustrates that there are both "red" 
and "blue" dusty star forming galaxies in clusters, and 
that our approach of modeling the 5.8/im and S.O/xm LFs 
with populations of both is reasonable. Furthermore, the 
fact that these are some of the brightest MIR sources in 
the cluster field, and that most have colors similar to 



the Sbc model, rather than the M82 model, is consis- 
tent with our conclusion that the 8.0/im LF at this red- 
shift is best modeled using the quiescent -l-regular model, 
with no need for a luminous dusty starburst component. 
We defer a more detailed discussion of the spectroscopy, 
including quantitative measurements of star formation 
from line widths to a future paper (Muzzin et al. 2008, 
in preparation). 

Overall, Figure 24 demonstrates that the galaxy tem- 
plates used to model the cluster LFs agree well with the 
colors of spectroscopically confirmed cluster galaxies, and 
that they are reasonable descriptions of star forming and 
non-star forming galaxies between < z < 1. 

6.4. Systematic Uncertainties 

The data presented in this paper support a self- 
consistent model of the evohition of stellar mass assem- 
bly and dusty star formation in clusters; however, there 
are several details of this analysis that have not been 
discussed and could potentially result in inappropriate 
conclusions being drawn from the data. Although it is 
difficult to quantify what effect, if any, these details will 
have on the interpretation of the data, we believe it is 
important to at least note these issues here. 

One worthwhile concern is the sample of clusters used 
in the analysis. Although this sample is much larger 
than the mere handful of clusters that have been studied 
in the MIR thus far, it is still of modest size and subject 
to cosmic variance. In particular, given that the clusters 
come from only 3.8 deg^, it is unclear whether the higher 
redshift clusters in the sample are truly the progenitors 
of the lower redshift clusters. Unfortunately, a cosmolog- 
ically significant sample of clusters covering of the order 
100 degree^ or more is likely needed to avoid biases that 
might result from cosmic variance in the sample. 

Another potential problem is that there are many more 
low richness clusters in the sample than high richness 
clusters, simply because of the nature of the cluster mass 
function. Any effects that depend on cluster mass will 
clearly be missed by combining these samples. This could 
be important because processes that could quench star 
formation (e.g., ram- pressure stripping, gas strangula- 
tion) or incite starbursts (tidal effects, harassment) will 
likely depend on cluster mass. Using a much larger sam- 
ple which can be separated by both mass and redshift 
would be invaluable for studying this issue further. 

Perhaps the most important concern is that there is 
a degeneracy between the intensity of star formation in 
clusters and the fraction of star forming galaxies. We 
showed in §6.3 that the color models used for the clus- 
ter galaxies reproduce the colors of cluster galaxies with 
spectroscopic redshifts very well; however, even though 
these colors are correct, the models of the 5.8/im and 
8.0/um LFs still depend on the assumed fb as a func;tion 
of magnitude and redshift for the clusters. If the f;, val- 
ues are overestimated and need to be reduced, then a 
larger fraction of dusty starburst galaxies than we have 
assumed will be required to correctly model the cluster 
5.8/im and 8.0/im LFs. Likewise, if the f;, is imderes- 
timated, fewer dusty starbursts will be required. The 
assumed f^ are consistent with most previous studies; 
however, optimally, if more data were available the fb 
should be calculated from the clusters themselves and 
this would avoid this degeneracy. 
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Fig. 24. — Plot of optical - IRAC or IRAC - IRAC colors of galaxies as a function of redshift. The red, purple, and green points are 
spectroscopic cluster members classified as non-star forming, star forming, and AGN respectively. The solid, dotted, and dashed lines are 
the model colors from the passive evolution model, the Sbc model, and the M82 model respectively. 



Lastly, it is worth mentioning that much of the excess 
seen in the 8.0/.tm LFs is near the hmiting magnitude of 
the survey. Problems with the background estimation 
could artificially inflate these values. It is unlikely that 
this is the case because if the excess of galaxies near the 
faint limit of the survey were due to an undersubtraction 
of the background, it should also be seen in the lower red- 
shift LFs, which it is not. Furthermore, undersubtraction 
of the background should be even more prevalent in the 
lower redshift LFs because clusters have much larger an- 
gular sizes and therefore more total area from which to 
undersubtract the background. It is unlikely that this is 
a problem; however, deeper data would be useful in en- 
suring there are no errors due to completeness near the 
survey limit. 

7. CONCLUSIONS 

We have presented a catalogue of 99 candidate clusters 
and groups at 0.1 < z^hot < 1-3 discovered in the Spitzer 
First Look Survey using the cluster red sequence tech- 
nique. Using spectroscopic redshifts from FLS foUowup 
campaigns and our own spectroscopic foUowup of clusters 
we have shown that the R - 3.6/im color of the cluster 
red sequence is an accurate photometric redshift estima- 
tor at the Az = 0.04 level at z < 1.0. Furthermore, we 
demonstrated that the properties of the FLS cluster cat- 
alogue are similar to previous cluster surveys such as the 
RCS-1. Using this cluster sample we studied the evolu- 
tion of the cluster 3.6/im, 4.5/im, 5.8/im, and 8.0/im LFs. 
The main results from these LFs can be summarized as 
follows: 



• In agreement with previous work, the evolution of 
the 3.6/im and 4.5/im LFs between 0.1 < z < 1.0 is 
consistent with a passively evolving population of 
galaxies formed in a single-burst at z > 1.5. Given 
that the 3.6/im and 4.5/im bandpasses are reason- 
able proxies for stellar mass, this suggests that the 
majority of stellar mass in clusters is already as- 
sembled into massive galaxies by 2 ~ 1. 

• The MIR color cuts used to select AGN by Lacy et 
al. (2004) and Stern et al. (2005) suggest that the 
fraction of cluster galaxies that host MIR-bright 
AGN at z < 0.7 is low. We estimate that the 
AGN fraction of cluster galaxies detected at 3.6/im 
is lli%. AGN are a larger, but still modest com- 
ponent of the 5. Sum and 8.0/im cluster population, 
approximately 3_3% of these galaxies. 

• The cluster 5.8/im and 8.0/im LFs do not look sim- 
ilar to the 3.6/im and 4.5/im LFs, and this is due 
to the presence of the cluster star forming galaxies. 
Star forming galaxies are much brighter in these 
bandpasses than early type galaxies and their vary- 
ing fractions with redshift cause deviations from 
the shape of the 3.6/im and 4.5/im LFs. The 5.8/im 
and 8.0/im LFs are well-described using different 
fractions of three basic types of galaxies: quies- 
cent, regular star forming, and dusty starburst by 
assuming that the fractions of the latter two are 
proportional to the cluster ff,. 

• The 8.0/im cluster LFs suggest that both the fre- 
quency and SSFR of star forming cluster galaxies 
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is increasing with increasing redshift. In particu- 
lar it appears that when compared to star forming 
galaxies in the local universe, the intensity of star 
formation in clusters evolves from "weak" to "regu- 
lar" to "star burst" with increasing redshift. Quali- 
tatively, this evolution mimics the evolution in the 
universal star formation density with redshift sug- 
gesting that this evolution is at least in part caused 
by the accretion of star forming galaxies into the 
cluster environment. 

• Comparing the 3.6y^m and 4.5y^m cluster and field 
LFs with similar normalization shows that the LFs 
in these environments arc similar, with evidence 
for a small excess in the brightest galaxies in clus- 
ters, likely caused by the cluster giant ellipticals. In 
agreement with previous K-band studies this sug- 
gests that the distribution of galaxies as a func- 
tion of stellar mass in both environments is roughly 
equivalent. 

• There is a significant differential evolution in the 
cluster and field 5.8/um and 8.0/um LFs with red- 
shift. At z = 0.65 the cluster is more abundant in 
8.0/xm galaxies than the field; however, thereafter 
the relative number of 5.8/im and S.Ofim galaxies 
declines in clusters with decreasing redshift and by 
^; = 0.15 the cluster is underdense in these sources 
by roughly a factor of 5. This differential evolution 
could be explained if starbursts are preferentially 
triggered during the early formation stages of the 
cluster but then preferentially quenched thereafter 
by the high density environment. 

A well-sampled spectroscopic study of several high- 

redshift clusters with MIR data would be extremely valu- 
able for verifying our interpretation of the IRAC cluster 
LFs because it is always difficult to draw incontrovert- 
ible conclusions from LFs alone. Still, the cluster LFs 
do show a strong increase in the number of 5.8/im and 
8.0/Lim sources in clusters with increasing redshift which 
must almost certainly be attributed to increased amounts 
of dusty star formation in higher redshift clusters. 

One of the strengths of this analysis is that it is based 
on a relatively large sample of galaxy clusters. It has 
become clear from the handful of clusters studied thus 
far by ISO and Spitzer that the MIR properties of clus- 
ter galaxies can be quite different from cluster to cluster. 
They may depend on dynamical state, mass, f(,, or other 



parameters (e.g., Coia et al. 2005; Gcach ct al. 2006). 
The advantage of using many clusters is that it provides 
a metric of how the "average" cluster is evolving as a 
function of redshift. Detailed studies of individual clus- 
ters with significant ancillary data will pave the way to a 
better understanding of the physics behind the evolution 
of dusty star formation in cluster galaxies; however, large 
statistical studies such as this one will indicate whether 
the clusters studied in future work are representative of 
the cluster population as a whole, or are potentially rare, 
biased clusters with unusual properties caused by an on- 
going merger or some other event. 

It is worth noting that although the quality of the LFs 
provided by the 99 clusters in the FLS is good, these LFs 
would still benefit from a larger statistical sample. In 
particular, a larger sample would allow for the separation 
of clusters by other properties such as mass or morphol- 
ogy, and to understand if these properties play a role in 
shaping the MIR cluster galaxy population. We are cur- 
rently working on a survey to detect clusters in the much 
larger SWIRE survey: the Spitzer Adaptation of the Red 
sequence Cluster Survey (SpARCS). This project has 13 
times more area than the FLS and is a factor of 2 deeper 
in integration time in the IRAC bands. The analysis of 
that sample should provide a significant improvement in 
the quality of the cluster LFs. 
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TABLE 1 — Continued 
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R.A. 


Decl. 






M200 


R-200 


Centroid 














J2000 


J2000 




M0 X 10^^ 


Mpc 




(1) 




(2) 




(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


1 o 


r ijo 


J1722094 


-5935.2 


U. I o 




1 7-99-nQ A 


oy .00. lu.o 


ouu 


J-UU 


1 on 

l.UU 


U.OU 


JTLo-r iux 


74 


FLS 


J1720354 


-5928.6 


0.78 




17'20:35.5 


^Q-^R-ACi A 
u ^? . iiO . L\/ . 1 


428 


176 


1.32 


0.94 


BCG 


75 


FLS 


J17141H 


-6027.7 


0.78 




17:14:11.7 


fin-27-.i 1 3 


705 


228 


2.97 


1.22 


BCG 


76 


FLS 


J1715454 


-5853.8 


0.78 




17:15:45.9 


■.JO . tj!) . 10 . \) 


291 


142 


0.71 


0.76 


BCG 


77 


FLS 


J1715564 


-5859.9 


n 7Q 




1 7. 1 c;. c;ft 1 
J- 1 . -LO. ou. -L 




uou 


91 7 


2 52 






78 




J1719324 


-5929.3 


n 7Q 




17-1 Q-'^O n 


oy .zy. ±0.0 


AQQ 
4^yy 


ly 1 


1 7n 

1. / U 


1 02 


rto-x iux 


79 


FLS 


J172019+5926.6 


0.79 




± 1 .zu. ±y.o 




291 


142 


0.71 


0.76 


J\,0 iT ILIA. 


80 


FLS 


J171828+5836.2 


0.79 




1 7'^R''?R 7 


uo.ou. j.0.0 


498 


191 


1.69 


1.02 




SI 

O -L 


FT S 


J172304+5832.3 


n SI 

U.O -L 






KQ.QO.I C « 
Oo.OZ. J.O.U 


ouo 


1 fii 


1 01 


U.OU 




82 


FLS 


J171657+6004.8 


0.82 




X 1 . xvj.u 1 . y 


fiO'flil'AQ ^ 
uu.u'i.^y.o 


711 


229 


3.02 


1.23 


BCG 


OO 


FT S 


J171945+5909.1 






1 7-1 Q'A^ A 


oy.uy.uy. -L 


ou 1 


1 0*^ 


1 7 A 


-L.UO 


xvo-ir lux 




FT S 


J171155+6013.1 


n Qo 

u.yu 




1 7-1 1 -f^^ 1 
1 / . ii .00. 1 


UU. J-O.Uo.O 


522 


1 Qfi 


1 82 


1 04 


xxo-r lux 




FLS 


,11718084 


-5915.8 


n Qi 




17-1 s-ns 7 


oy . J-O.ovj. ( 


00 1 


1 fiS 


1 12 


u.oy 


rvo-f iux 


OO 


FLS 


,11712234 


-6015.1 


n 

u. yo 




1 7.1 9-9Q 7 


fin- 1 f^-nQ A 


ouo 


211 


2 32 


1 1 


R*:i PliiY 
rvo-x iiix 


87 


FLS 


J17105H 


-5930.8 


1.02 




i7-in-'^i R 

L 1 . l.\J .0± .0 


'^Q '^O '^n 
oy .ou . ou . 


760 


237 


3.36 


1.27 


BCG 


88 


FLS 


J1721474 


-6011.5 


1.02 




17:21:47.3 


uu. X ± . 00. 1 


277 


141 


0.65 


0.74 


BCG 


89 


FLS 


J1718524 


-6009.9 


1.02 




17:18:52.7 


60:09:56.9 


485 


189 


1.62 


1.00 


RS-Flux 


90 


FLS 


J17122H 


-6010.6 


1.03 




17:12:21.2 


60:10:41.0 


349 


160 


0.95 


0.84 


BCG 


91 


FLS 


J17143H 


-5946.9 


1.06 




17:14:31.9 


59:46:59.5 


425 


177 


1.31 


0.94 


RS-Flux 


92 


FLS 


J1711174 


-5902.8 


1.06 




17:11:17.5 


59:02:48.6 


287 


144 


0.69 


0.76 


BCG 


93 


FLS 


J1721264 


-5856.6 


1.11 


1.045(4) 


17:21:26.4 


58:56:41.7 


646 


218 


2.58 


1.17 


BCG 


94 


FLS 


J1712274 


-6015.2 


1.14 




17:12:27.0 


60:15:16.7 


448 


182 


1.43 


0.96 


BCG 


95 


FLS 


J1720454 


-5834.8 


1.17 




17:20:45.3 


58:34:50.9 


386 


169 


1.12 


0.89 


RS-Flux 


96 


FLS 


J172113+5901.0 


1.24 




17:21:13.8 


59:01:05.7 


338 


158 


0.90 


0.83 


RS-Flux 


97 


FLS 


J171223+6006.9 


1.27 




17:12:23.6 


60:06:56.4 


208 


124 


0.41 


0.64 


RS-Flux 


98 


FLS 


J171942+5938.3 


1.38 




17:19:42.8 


59:38:23.2 


374 


165 


1.06 


0.87 


RS-Flux 



Note. — (3) Photometric redshift estimated from red sequence color, (4) Mean spectroscopic redshift of galaxies with red sequence 
weights > 0.2, the number of spectroscopic rcdshifts is included in brackets, (7) Cluster richness parameterized by Bgc,i{, (8) Error in 
Bgc,i?> (9) M200 estimated from Bgc,H using eq. 1. (10) R200 estimated from Bg^^ij using eq. 2., (11) Best centroid of the cluster. 



TABLE 2 
Assumed Blue Fractions 



z 


F), (M < M*) 


(1) 


(2) 


0.15 
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